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Abstract 

Segmentation  and  crustal  accretion  at  oceanic  transform  fault  systems  are 
investigated  through  a  combination  of  geophysical  data  analysis  and  geodynamical  and 
geochemical  modeling.  Chapter  1  examines  the  effect  of  fault  segmentation  on  the 
maximum  predicted  earthquake  magnitude  of  an  oceanic  transform  fault  system.  Results 
of  thermal  modeling  suggest  that  fault  segmentation  by  intra-transform  spreading  centers 
(ITSC)  drastically  reduces  the  available  brittle  area  of  a  transform  fault  and  thus  limits 
the  available  earthquake  rupture  area.  Coulomb  stress  models  suggest  that  long  ITSCs 
will  prohibit  static  stress  interaction  between  segments  of  a  transform  system  and  further 
limit  the  maximum  possible  magnitude  of  a  given  transform  fault  earthquake.  In  Chapter 
2,  gravity  anomalies  from  a  global  set  of  oceanic  transform  fault  systems  are  investigated. 
Surprisingly,  negative  residual  mantle  Bouguer  gravity  anomalies  are  found  within  fast¬ 
slipping  transform  fault  domains.  These  gravity  observations  suggest  a  mass  deficit 
within  fast-slipping  transform  faults,  which  may  result  from  porosity  variations,  mantle 
serpentinization,  and/or  crustal  thickness  variations.  Two-dimensional  forward  modeling 
and  the  correlation  of  the  negative  gravity  anomalies  to  bathymetric  highs  indicate  crustal 
thickness  excesses  in  these  locations.  Finally,  in  Chapter  3,  mantle  thermal  and  melting 
models  for  a  visco-plastic  rheology  are  developed  to  investigate  the  process  of  mantle 
melting  and  crustal  accretion  at  ITSCs  within  segmented  transform  faults,  and  are  applied 
to  the  Siqueiros  transform  fault  system.  Models  in  which  melt  migrates  into  the 
transform  fault  domain  from  a  large  region  of  the  mantle  best  explain  the  gravity-derived 
crustal  thickness  variations  observed  at  the  Siqueiros  transform.  Furthermore,  a  mantle 
potential  temperature  of  1350°C  and  fractional  crystallization  at  depths  of  9  -  15.5  km 
best  explain  the  major  element  composition  variation  observed  at  the  Siqueiros  transform. 
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Introduction 


The  global  mid-ocean  ridges  are  segmented  by  more  than  120  transform  fault 
systems.  These  transform  fault  systems  are  critical  for  accommodating  the  geometry  of 
the  ridge  system  as  well  as  crustal  accretion.  The  primary  focus  of  this  Ph.D.  thesis  is 
the  investigation  of  oceanic  transform  faults  using  a  multidisciplinary  approach.  Chapter 
1  explores  the  seismicity  of  oceanic  transform  faults  and  examines  the  relationship 
between  fault  segmentation  and  earthquake  rupture.  Chapter  2  utilizes  residual  mantle 
Bouguer  gravity  anomaly  (RMBA)  calculations,  3D  thermal  models,  and  2D  forward 
models  to  investigate  the  structure  of  a  global  set  of  transform  faults.  Finally,  in  Chapter 
3,  1  develop  combined  geodynamical  and  geochemical  models  to  constrain  melt  generation 
and  extraction  at  segmented  transform  fault  systems.  Below,  1  present  a  brief  description 
of  each  of  these  studies  and  their  primary  conclusions. 

In  Chapter  1,  I  combine  Coulomb  static  stress  models  with  available  geological  and 
seismicity  data  to  illustrate  the  importance  of  fault  segmentation  in  models  of  earthquake 
interaction  at  oceanic  transform  faults.  The  2D  thermal  models  demonstrate  that  fault 
segmentation  decreases  the  maximum  magnitude  of  earthquakes  by  decreasing  the  brittle 
area  thus  limiting  the  fault  rupture  of  the  transform  fault  system  [ Gregg  et  al. ,  2006], 
Coulomb  static  stress  models  are  used  to  explore  whether  adjacent  fault  segments  behave 
independently  of  each  other,  and  how  their  stress  interaction  depends  on  their  offset 
distance.  I  find  that,  if  intra-transform  spreading  centers  are  sufficiently  long,  adjacent 
transform  fault  segments  will  be  decoupled  and  behave  independently  of  each  other.  I 
also  demonstrate  the  likelihood  of  stress  interaction  between  strike-slip  earthquakes  and 
normal  faults  at  nearby  spreading  centers.  These  findings  are  important  for  understanding 
the  seismic  budget  of  oceanic  transform  faults. 

In  Chapter  2,  I  analyze  3D  gravity  anomalies  for  19  oceanic  transform  fault 
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systems  globally.  By  evaluating  transform  faults  from  several  mid-ocean  ridge  systems 
using  3D  gravity  inversions,  3D  thermal  models,  and  2D  forward  models,  1  show  that 
there  is  systematic  variation  in  the  gravity  signature  of  oceanic  transform  faults  as  a 
function  of  their  slip  rate  [Gregg  et  al.,  2007b].  The  spreading-rate  dependency  of  the 
gravity  anomalies  is  especially  evident  at  fast-slipping  transform  faults  where  I  find  more 
negative  RMBA  values  in  comparison  to  the  adjacent  ridge  axes.  All  previous  studies  of 
slow-spreading  mid-ocean  ridges  have  shown  more  positive  RMBA  values  within  the 
transform  fault  domain  instead,  indicating  crustal  thinning  at  the  end  of  ridge  segments  and 
focused  crustal  accretion  at  ridge  segment  centers  [Kuo  and  Forsyth ,  1988;  Lin  et  al., 
1990],  Based  on  2D  forward  models  and  available  geological  data,  I  propose  that  the 
negative  RMBA  values  observed  at  fast-slipping  transform  faults  indicate  that  fast¬ 
slipping  systems  have  excess  magma  supply  and  increasing  crustal  production  within  the 
transform  fault  domain.  These  findings  directly  contradict  the  classic  models  of  crustal 
accretion  at  mid-ocean  ridges  and  require  further  investigation  to  determine  the  source  of 
the  negative  RMBA  at  fast-slipping  transform  faults. 

In  the  final  chapter  of  my  thesis,  I  look  at  the  process  of  mantle  melting  and  melt 
extraction  at  fast-slipping  transform  faults  in  an  effort  to  explain  the  findings  of  increased 
crustal  accretion  at  these  systems  [Gregg  et  al.,  2007b],  I  combine  thermal  models  of 
segmented  ridge-transform-ridge  systems  with  fractional  melting  [Kinzler  and  Grove, 
1992b;  a;  1993]  and  fractional  crystallization  models  [Yang  et  al.,  1996].  The  melt 
migration  model  of  Sparks  and  Parmentier  [1991]  is  used  to  investigate  the  source  region 
of  intra-transform  spreading  centers  and  examine  two  end-member  melt  aggregation 
models  of  wide  vs.  narrow  melt  pooling  regions.  I  apply  these  models  to  the  specific  case 
of  Siqueiros  transform  fault  on  the  East  Pacific  Rise  where  there  has  been  extensive 
geological  and  geophysical  mapping  [Fornari  et  al.,  1989;  Carbotte  and  Macdonald,  1992; 
Gregg  et  al.,  2007b],  as  well  as  sampling  via  dredge  and  Alvin  dives  [Perfit  et  a!.,  1996; 
Flays  et  al.,  2004],  Our  findings  suggest  that  in  order  to  explain  the  estimated  crustal 
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thicknesses  indicated  by  the  RMBA  calculated  at  Siqueiros,  melt  must  be  pooled  from 
large  distances  into  the  transform  fault  domain.  Furthermore,  we  find  that  a  mantle 
potential  temperature  of  1350C  and  fractional  crystallization  at  depths  of  9-15.5  km  best 
fits  the  geochemical  observations  at  the  Siqueiros  transform  fault  system  [Perfit  et  al. , 
1996;  Hays,  2004].  These  melt  migration  models  provide  us  with  first-order  insights  into 
the  geological  processes  that  may  be  occurring  at  the  Siqueiros  transform  fault. 

Future  Directions 

The  research  contained  in  this  thesis  is  a  starting  point  for  the  development  of 
many  exciting  projects  for  the  coming  years.  In  the  future,  I  plan  to  extend  my 
investigation  to  study  several  first-order  geological  processes. 

Coulomb  stress  modeling  is  an  invaluable  tool  for  estimating  static  stress  transfer 
following  an  earthquake  or  diking  episode  and  determining  likely  locations  for  subsequent 
events  [e.g.,  Kinget  at. ,  1994],  As  part  of  this  thesis,  I  successfully  used  stress  models  to 
illustrate  earthquake  interaction  at  oceanic  transform  faults  and  the  stress  interaction 
between  strike-slip  earthquakes  and  normal  faults  at  nearby  spreading  centers  [ Gregg  et 
a/.,  2006].  It  has  also  been  suggested  that  stress  transfer  due  to  tectonic  events  can 
stimulate  interactions  between  volcanoes  and  earthquakes.  In  particular,  recent  research 
on  Hawaii  has  illustrated  striking  correlations  between  dike  intrusions  and  earthquake 
occurrence  at  Mauna  Loa  [Walter  and  Amelnng,  2004;  2006;  Amelung  et  al.,  2007; 
McGovern,  2007],  These  techniques  can  also  be  applied  to  understanding  continental 
rifting;  specifically,  what  is  the  stress  interaction  between  dike  injection  during  rifting 
change  and  subsequent  faulting  events?  For  example,  temporal  and  spatial  seismicity 
patterns  along  the  East  African  Rift  system  indicate  that  there  is  ongoing  stress 
interaction  between  adjacent  rift  segments  via  diking  and  faulting  events  [Aronovitz  et  al ., 
2007],  Increasing  our  understanding  of  how  earthquakes  and  volcanoes  interact  in  settings 
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such  as  Hawaii  and  Africa  is  also  important  for  aiding  our  ability  to  assess  future  hazards, 
which  is  a  direction  1  plan  to  pursue. 

The  spreading  rate  dependence  of  mid-ocean  ridge  morphology  is  widely 
observed.  Typical  slow-spreading  ridges  are  characterized  by  axial  valleys  and  large 
topographic  variations,  whereas  fast-spreading  ridges  are  generally  associated  with  axial 
highs  and  smoother  topography  [ Macdonald ,  1983],  However,  in  many  locations  it 
appears  that  magma  supply  may  play  an  even  greater  role  in  controlling  the  ridge 
morphology  than  spreading  rate.  This  paradox  remains  one  of  the  first-order  problems  in 
understanding  the  dynamics  of  mid-ocean  ridges.  For  example,  the  Cleft  segment  of  the 
Juan  de  Fuca  Ridge  north  of  the  350  km-long  Blanco  transform  fault  has  the  morphology 
of  a  fast-spreading  ridge  and  is  thought  to  be  magmatically  robust,  while  the  Gorda  ridge 
to  the  south  of  Blanco,  which  is  also  spreading  at  a  rate  of  5.8  cm/yr,  exhibits  an  axial 
valley  typical  of  a  slow-spreading  ridge  and  is  considered  magmatically  starved  [Embley 
and  Wilson ,  1992;  Hooft  and  Detrick,  1995],  Other  examples  include  two  adjacent 
segments  at  ~18°S  on  the  fast-spreading  East  Pacific  Rise  separated  by  a  small  non¬ 
transform  discontinuity  [Sinton  el  al,  2002]  and  adjacent  segments  at  ~1 14°E  along  the 
intermediate  spreading  South  East  Indian  Ridge  (SEIR),  which  are  separated  by  a  100 
km-long  transform  fault  [Cochran  and  Sempere,  1997;  Buck  et  al.,  2005;  Gregg  et  al., 
2007b],  Not  only  is  the  apparent  decrease  in  magma  budget  visible  in  the  morphology  of 
these  ridge  segments,  it  is  also  observed  in  seismic  and  gravity  data  [Hooft  and  Detrick, 
1995]  as  well  as  in  the  major  and  trace  element  geochemistry  of  lavas  from  these  settings 
[Rubin  and  Sinton,  2007],  At  present  it  is  unclear  why  adjacent  segments  spreading  at  the 
same  rate  are  in  such  stark  contrast  to  each  other,  and  our  current  understanding  of  the 
physics  behind  mid-ocean  ridge  spreading  does  not  fully  address  this  problem.  Possible 
explanations  may  include  thermal  effects  due  to  ridge  offsets,  regional  tectonics,  mantle 
flow  variations,  mantle  buoyancy  [Buck  et  al.,  2005],  and  mantle  source  heterogeneities.  I 
plan  to  explore  this  research  direction  in  the  future. 
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In  chapter  3  of  this  thesis  I  have  developed  techniques  for  coupling  geodynamical 
and  geochemical  models,  which  can  be  constrained  by  geochemical  and  geophysical  data  to 
examine  the  driving  forces  governing  the  magma  budget  of  a  ridge  segment  and  the 
resultant  ridge  morphology  and  geochemistry.  These  models  can  be  further  developed  by 
accounting  for  the  effect  of  fractional  melting  and  crystallization  on  the  trace  element 
composition  of  extruded  lavas  at  mid-ocean  ridges.  As  part  of  my  thesis  work,  we  have 
made  significant  advancements  in  understanding  geodynamical  models  by  coupling  them 
with  fractional  melting  [Kinzler  and  Grove,  1992a,  b;  1993]  and  fractional  crystallization 
models  [Yang  et  a!.,  1996].  Incorporation  of  trace  elements  into  these  models  will  further 
strengthen  their  capabilities  and  provide  another  independent  variable  for  comparison 
with  geochemical  observations,  a  direction  that  I  am  interested  in  exploring  in 
collaboration  with  geochemist  colleagues. 

Hydrated  oceanic  lithosphere  is  an  important  component  of  the  subduction  zone 
system.  It  is  widely  accepted  that  water  plays  a  fundamental  role  in  the  processes  of 
mantle  melting  and,  subsequently,  arc  volcanism  in  the  subduction  setting.  However, 
very  little  is  known  about  the  actual  water  budget  available  from  hydration  of  the  oceanic 
lithosphere.  Two  popular  models  for  oceanic  lithosphere  hydration  include:  1)  the 
serpentinization  of  the  upper  lithosphere  at  oceanic  transform  faults  and  fracture  zones; 
and  2)  hydration  of  the  upper  mantle  to  depths  of  40  km  at  large-scale  normal  faults 
flanking  the  outer  rise  (thought  to  be  the  result  of  the  bending  of  the  slab  as  it  subducts) 
[Hacker  et  al.,  2003].  However,  it  is  unclear  whether  these  two  mechanisms  are  viable 
and  if  so,  how  much  water  is  locked  into  hydrous  phases  within  the  down-going  plate  at 
these  locations.  This  is  another  direction  for  future  exploration. 

During  my  thesis  work,  I  have  become  increasingly  interested  in  the  effects  of 
lithospheric  hydration  as  a  way  to  explain  negative  gravity  anomalies  found  within  fast¬ 
slipping  oceanic  transform  faults  [Gregg  et  al.,  2007b]  as  well  as  seismic  moment  deficits 
observed  at  many  oceanic  transform  faults  [Boettcher  and  Jordan ,  2004],  In  the  past 
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couple  of  years  we  have  developed  several  3D  thermal  models  for  various  transform  fault 
geometries  and  mantle  rheologies  [Behn  et  al.,  2007;  Gregg  et  al.,  2007a],  as  well  as  the 
inclusion  of  hydrous  phase  stability  fields  to  assess  localized  hydration  [ Roland  et  al., 
2007].  A  powerful  result  of  these  modeling  efforts  is  the  ability  to  estimate  water 
budgets  expected  from  oceanic  transform  faults  of  varying  spreading  rate  and  geometry. 
The  next  step  in  these  calculations  is  to  look  at  hydration  in  the  inactive  fracture  zones 
and  explore  the  dehydration  of  these  sections  of  the  oceanic  lithosphere  as  they  are 
incorporated  at  subduction  zones  at  settings  such  as  the  Chile  Triple  Junction  or  the 
Blanco  Fracture  Zone. 
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Chapter  1:  Segmentation  of  transform  systems  on  the  East  Pacific  Rise: 
Implications  for  earthquake  processes  at  fast-slipping  oceanic  transform 
faults 


Abstract 

Seven  of  the  eight  transform  systems  along  the  equatorial  East  Pacific  Rise  from 
12°N  to  15°S  have  undergone  extension  due  to  reorientation  of  plate  motions  and  have 
been  segmented  into  two  or  more  strike-slip  fault  strands  offset  by  intratransform 
spreading  centers  (ITSCs).  Earthquakes  recorded  along  these  transform  systems  both 
teleseismically  and  hydroacoustically  suggest  that  segmentation  geometry  plays  an 
important  role  in  how  slip  is  accommodated  at  oceanic  transforms.  Results  of  thermal 
calculations  suggest  that  the  thickness  of  the  brittle  layer  of  a  segmented  transform  fault 
could  be  significantly  reduced  by  the  thermal  effect  of  ITSCs.  Consequently,  the  potential 
rupture  area,  and  thus  maximum  seismic  moment,  is  decreased.  Using  Coulomb  static 
stress  models,  we  illustrate  that  long  ITSCs  will  prohibit  static  stress  interaction  between 
transform  segments  and  limit  the  maximum  possible  magnitude  of  earthquakes  on  a  given 
transform  system.  Furthermore,  transform  earthquakes  may  have  the  potential  to  trigger 
seismicity  on  normal  faults  flanking  ITSCs. 


Published  as  P.  M.  Gregg,  J.  Lin,  and  D.  K.  Smith,  Segmentation  of  transform  systems  on  the  East  Pacific 
Rise:  Implications  for  earthquake  processes  at  fast-slipping  oceanic  transform  faults.  Geology,  v.  34,  no.  4. 
p.  289-292,  2006. 
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Introduction 

Segmented  transform  systems  are  composed  of  several  fault  strands  offset  by  short 
ridges  or  rifts  referred  to  as  intra-transform  spreading  centers  (ITSCs)  (Menard  and 
Atwater,  1969;  Searle,  1983;  Pockalny  et  al.,  1997),  where  active  seafloor  spreading  and 
crustal  accretion  are  occurring  (Fomari  et  al.,  1989;  Hekinian  et  al.,  1992;  Perfit  et  al., 
1996).  Along  the  equatorial  East  Pacific  Rise  between  15°S  and  12°N  (Fig.  1),  the 
Siqueiros,  Quebrada,  Discovery,  Gofar,  Yaquina,  Wilkes,  and  Garrett  transform  systems 
have  all  undergone  transtension  due  to  changes  in  plate  motions,  and  each  of  these 
transforms  is  segmented  by  at  least  one  1TSC  (Searle,  1983;  Fomari  et  al.,  1989; 
Lonsdale,  1989;  Goff  et  al.,  1993;  Pockalny  et  al.,  1997).  The  Clipperton  transform 
system  has  undergone  several  periods  of  transpression  (Pockalny,  1997),  and  is  the  only 
unsegmented  transform  system  along  the  equatorial  East  Pacific  Rise. 

The  global  deficiency  of  seismic  moment  release  on  oceanic  transform  systems  has 
led  researchers  to  hypothesize  that  a  significant  portion  of  oceanic  transform  slip  is 
accommodated  aseismically  (e.g.,  Boettcher  and  Jordan,  2004).  However,  global 
seismicity  studies  have  yet  to  consider  the  prevalence  of  transform  fault  segmentation. 
Dziak  et  al.  (1991)  observed  that  earthquake  sizes  generally  correlate  with  the  lengths  of 
individual  fault  segments  at  the  Blanco  transform  fault.  Our  observations  of  earthquakes 
recorded  on  East  Pacific  Rise  transform  faults  indicate  that  segmentation  is  an  important 
factor  influencing  rupture  of  large  earthquakes  at  oceanic  transforms.  While  it  has  been 
shown  that  segmentation  and  fault  steps  play  an  important  role  in  controlling  the 
earthquake  behavior  of  continental  strike-slip  faults  (e.g.,  Harris  and  Day,  1993),  the 
influence  of  segmentation  and  ITSCs  on  earthquake  processes  at  an  oceanic  transform 
system  has  not  been  studied  in  detail. 

In  this  paper,  we  use  teleseismically  and  hydroacoustically  recorded  seismicity  data 
from  the  equatorial  East  Pacific  Rise  and  Coulomb  static  stress  models  to  explore  the 
effect  of  ITSCs  on  static  stress  interaction  between  transform  fault  segments.  We 
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Figure  1.  Regional  map  of  the  equatorial  EPR  showing  large  transform  and  non¬ 
transform  offsets.  Segmentation  geometry  is  included  based  on  previous  geological 
mapping  of  the  transform  systems  (e.g.,  Fomari  et  al.,  1989;  Lonsdale,  1989).  Inset: 
Regional  map  showing  location  of  the  full  array  of  NOAA  Pacific  Marine  Environmental 
Laboratory  hydrophones. 
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Figure  2.  Comparison  of  estimated  areas  of  brittle  lithosphere  using  a  one-dimensional, 
steady-state  lithosphere  cooling  model  (McKenzie,  1969)  for  the  Clipperton  (B)  and 
Siqueiros  (C)  transforms.  A:  The  90  km  Clipperton  transform  system  (X-X’)  and  the  150 
km  Siqueiros  transform  system  (Y-Y’),  which  is  broken  into  five  major  segments  SI,  S2, 
S3,  S4,  and  S5  separated  by  four  ITSCs  SA,  SB,  SC,  and  SD  (Fomari  et  al.,  1989).  B: 
Calculated  area  of  brittle  lithosphere  for  temperatures  <600  8C  (shaded  region)  for  the 
Clipperton  transform.  C:  Comparison  of  the  calculated  areas  of  brittle  lithosphere  for  the 
Siqueiros  transform  for  a  model  of  unsegmented  geometry  (area  above  the  dotted  line) 
versus  a  model  consisting  of  five  individual  segments  offset  by  steady-state  ITSCs 
(shaded  area). 
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investigate  whether  adjacent  fault  segments  can  behave  independently  of  one  another, 
and  how  the  interaction  between  segments  depends  on  their  offset  distance. 

Transform  Segmentation 

Segmentation  of  the  transtensional  transform  systems  at  the  equatorial  East  Pacific 
Rise  has  resulted  in  individual  strike-slip  fault  strands  with  lengths  of  18-89  km,  with  an 
average  of  ~37  km.  The  ITSCs  separating  the  fault  strands  have  lengths  of  5-20  km,  with 
an  average  length  of  ~  1 1  km.  Fresh  lavas  collected  from  the  ITSCs  within  the  Siqueiros 
and  Garrett  transforms  (Hekinian  et  al.,  1992;  Perfit  et  al.,  1996)  indicate  that  ITSCs  are 
magmatically  active,  implying  that  the  regions  beneath  them  are  hotter,  and  thus  the 
lithospheric  plate  is  thinner  than  the  surrounding  domains.  To  explore  the  effect  of 
segmentation  on  the  transform  fault  thermal  structure,  we  use  a  half-space  steady-state 
lithospheric  cooling  model  (McKenzie,  1969;  Abercrombie  and  Ekstrom,  2001).  The 
temperature  within  the  crust  and  mantle,  T,  is  defined  as  T=  7’,„erf[y(2K/)'  ],  where  Tm  is 
the  mantle  temperature  at  depth,  assumed  to  be  1300°C;  v  is  the  depth;  k  is  the  thermal 
diffusivity,  assumed  to  be  1,026  mV;  and  t  is  the  age  of  the  lithosphere  obtained  by 
dividing  distance  from  the  ridge  axis  by  half  the  spreading  rate. 

Figure  2  compares  the  calculated  areas  of  brittle  deformation,  defined  as  regions  with 
calculated  temperatures  <  600°C,  for  the  geometries  of  the  Clipperton  and  Siqueiros 
transform  systems.  The  calculated  area  under  the  600  °C  isotherm  for  the  Clipperton 
transform  fault  is  326  knT,  compared  to  710  knT  for  a  model  of  a  single  unsegmented 
fault  with  the  cumulative  length  of  the  Siqueiros  transform  system.  When  the  actual 
segmentation  geometry  of  the  Siqueiros  transform  system  is  considered,  however,  the 
integrated  area  of  the  calculated  brittle  deformation  region  is  decreased  by  -60%  to  277 
km2. 

Seismic  moment  (Mc),  which  reflects  the  energy  released  by  an  earthquake,  is  a 
function  of  the  rupture  area  of  the  fault.  Specifically,  M„  =  G  x  D  x  S,  where  G  is  the 
shear  modulus,  estimated  to  be  27  GPa  from  seismic  velocities  (Canales  et  al.,  2003),  D 
is  the  average  slip,  and  S  is  the  estimated  brittle  area.  The  resulting  moment  magnitude  is 
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Mw  =  (2/3)  x  log(M0)  -10.73.  For  a  model  of  constant  stress  drop  during  pure  strike-slip 
earthquakes,  M0  =  (7t/2)  xAoxh'xS,  where  Ao  is  the  earthquake  stress  drop,  and  w  is 
the  fault  width,  estimated  from  S  divided  by  fault  length  (Stein  and  Wysession,  2003). 

Curves  in  Figure  3  show  the  predicted  earthquake  magnitudes  for  a  given  fault  area, 
assuming  models  of  constant  fault  slip  (Fig.  3A)  or  constant  stress  drop  (Fig.  3B)  during 
earthquakes.  Earthquakes  recorded  teleseismically  as  listed  in  the  Harvard  Centroid- 
Moment  Tensor  (CMT)  catalog  were  relocated  to  a  specific  transform  segment  using 
hydroacoustically  recorded  earthquakes,  which  have  smaller  location  errors  (<6  km)  (Fox 
et  al.,  2001).  The  maximum  earthquake  magnitudes  observed  on  each  of  the  transform 
fault  segments  at  the  equatorial  East  Pacific  Rise  from  1996  to  2001  are  plotted  in  Figure 
3.  Assuming  the  complete  rupture  of  a  given  individual  fault  segment,  we  can  estimate 
the  amount  of  slip  or  the  stress  drop  for  a  given  earthquake.  For  example,  the  largest 
earthquake  observed  on  the  Clipperton  transform  has  a  Mw  of  6.6  (Fig.  3).  If  the  entire 
brittle  area  of  Clipperton  (326  km2)  was  ruptured  during  this  earthquake,  the  estimated 
average  slip  is  1  m,  or  the  estimated  average  stress  drop  is  53  bar  or  5.3  MPa. 

Coulomb  Stress  Calculations 

Evidence  for  potential  earthquake  interactions  along  oceanic  transform  faults  has 
been  noted  in  several  investigations  (e.g.,  Toda  and  Stein,  2000;  Bohnenstiehl  et  al., 
2002,  2004;  McGuire  et  al.,  2002;  Forsyth  et  al.,  2003).  We  utilize  the  methods  of  King 
et  al.  (1994)  and  Toda  et  al.  (1998)  to  calculate  how  static  stress  from  a  moderate-sized 
earthquake  is  transferred  to  adjacent  faults,  and  assess  the  likelihood  of  rupturing 
multiple  transform  segments  during  a  single  earthquake.  According  to  Coulomb  failure 
criteria,  when  an  earthquake  occurs  on  a  source  fault,  changes  in  Coulomb  failure  stress 
(AOf)  on  a  receiver  fault  are  expressed  as  AOf  =  Axs  +  p'  x  Aon,  where  Axs  and  Aon  are 
changes  in  shear  and  normal  stresses,  on  the  receiver  fault,  and  p'  is  the  apparent  friction 
coefficient  adjusted  for  the  pore  pressure  effect  (King  et  al.,  1994). 
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Figure  3.  A:  The  predicted  maximum  moment  magnitude,  Mw,  of  earthquakes  (curved 
lines)  for  a  given  transform  segment  area  and  a  constant  slip  of  0.1,  0.3,  and  1.0  m.  The 
calculations  assume  that  the  earthquake  ruptures  the  entire  transform  segment.  Black 
dots  mark  the  observed  maximum  Mw  recorded  on  the  transform  segments  of  the 
equatorial  East  Pacific  Rise.  The  rightmost  data  point,  Mw  =  6.6,  corresponds  to  the 
Clipperton  transform.  B:  The  predicted  maximum  Mw  assuming  a  constant  stress  drop  of 
transform  earthquakes  of  1,  10,  and  100  bar. 
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Figure  4.  Schematic  models  showing  the  geometry  of  two  transform  segments  bisected 
by  a  single  ITSC  of  variable  length,  L.  The  source  earthquake  is  located  on  the  bottom 
right  transform  segment  with  its  left  edge  located  at  a  distance,  d,  from  the  ITSC- 
transform  intersection.  The  source  earthquake  is  assumed  to  be  a  strike-slip  event  on  a 
vertical  plane  parallel  to  the  transform  segment.  A:  A  scenario  in  which  the  receiver  fault 
is  a  strike-slip  fault  located  on  the  top  left  transform  segment,  which  is  assumed  to  have 
the  same  dip,  strike,  and  rake  as  the  source  earthquake.  B:  A  scenario  in  which  the 
receiver  fault  is  a  normal  fault  located  along  the  ITSC,  which  is  assumed  to  have  a  dip  of 
608  and  is  parallel  to  the  ITSC. 
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We  consider  a  simple  geometry  in  which  two  adjacent  transform  segments  are  offset 
by  an  ITSC  of  variable  length,  L,  for  two  scenarios  assuming  the  receiver  faults  are  either 
strike-slip  events  along  the  adjacent  transform  segment  (Fig.  4A)  or  normal-faulting 
events  located  along  the  ITSC  (Fig.  4B).  The  rupture  size  for  the  source  earthquake  in 
both  cases  is  varied  to  reflect  typical  earthquake  magnitudes  observed  along  the 
segmented  transforms  of  the  equatorial  East  Pacific  Rise. 

For  the  first  scenario,  the  calculated  maximum  change  in  static  stress,  1 1  Aut  1 1 , 
transferred  to  the  receiver  fault  is  plotted  versus  L  for  Mw  =  5.0,  5.5,  and  6.0  (Fig.  5).  As 
the  separation  distance  between  the  two  transform  segments  increases,  the  predicted 
maximum  induced  Coulomb  stress  change  on  the  receiver  fault  decreases.  For  example,  if 
Mw  =  6.0,  d  =  0  km,  and  L  is  increased  from  5  to  15  km,  the  calculated  1 1  A<Jt  1 1  decreases 
from  1.35  bar  to  0.25  bar  (Fig.  5A).  The  proximity  of  the  earthquake  to  the  ITSC- 
transform  intersection  (1TI)  is  also  very  important:  the  closer  the  source  earthquake  is 
located  to  the  1T1  (i.e.,  smaller  d),  the  greater  the  predicted  maximum  Coulomb  stress 
change  on  the  receiver  fault  (Fig.  5).  For  example,  if  Mw  =  6.0,  L  =  5  km,  and  d  is 
increased  from  0  to  10  km,  the  calculated  ||Aof)|  decreases  from  2.5  bars  to  0.1  bars. 
Previous  studies  have  shown  statistically  significant  correlations  between  regions  of 
seismicity  rate  changes  following  continental  strike-slip  earthquakes  and  areas  of 
calculated  stress  changes  >  0.2- 1 .0  bars  (Toda  et  al.,  1 998). 

In  the  second  scenario  (Fig.  6),  for  Mw  >  5.0,  a  source  earthquake  with  relatively 
small  d  is  calculated  to  cause  significant  Coulomb  static  stress  increases  on  ITSC  parallel 
secondary  normal  faults.  Such  Coulomb  stress  changes  correspond  to  a  decrease  in 
normal  confining  pressure  across  the  ITSC  axis,  which  may  result  in  triggering  of  normal 
faulting  earthquakes  or  magmatic  diking  events  along  the  ITSCs.  The  predicted 
Coulomb  stress  changes  on  secondary  normal  faults  along  the  ITSC  are  a  strong  function 
of  the  location  of  the  source  earthquake.  For  example,  if  Mw  =  6.0  and  d  is  increased 
from  5  to  10  km,  the  calculated  ||Acf||  decreases  from  2.5  bars  to  0.5  bar. 
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Figure  5.  Calculated  maximum  Coulomb  stress  changes  on  a  secondary  strike-slip 
receiver  fault  caused  by  a  strike-slip  source  earthquake  (geometry  shown  in  Fig.  4A)  for 
source  earthquake  Mw  5  6.0,  5.5,  and  5.0.  Note  that  the  vertical  scale  is  different  for  each 
panel.  All  stress  calculations  were  carried  out  using  a  three-dimensional  boundary- 
element  model,  Coulomb  2.6  (Toda  et  al.,  1998),  assuming  that  both  the  source  and 
receiver  faults  extend  to  a  depth  of  5  km.  For  each  Coulomb  calculation,  we  used  a 
Young’s  modulus  of  62.5  GPa,  a  Poisson’s  ratio  of  0.25,  an  apparent  friction  coefficient, 
m9,  of  0.4  (e.g.,  King  et  ah,  1994),  and  a  tapered  slip  distribution.  Stresses  are  sampled 
on  a  horizontal  plane  at  a  depth  of  2  km.  The  maximum  change  in  Coulomb  stress  is 
taken  directly  from  the  point  on  the  receiver  fault  where  Coulomb  stress  reaches  a 
maximum  value.  Calculations  were  carried  out  for  L  5  1  to  20  km  and  d  5  0,  2.5,  5,  and 
10  km.  ITSC. 
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Figure  6.  Calculated  maximum  Coulomb  stress  changes  on  a  secondary  normal  receiver 
fault  along  the  ITSC  caused  by  a  strike-slip  source  earthquake  on  the  adjacent  transform 
segment  with  geometry  shown  in  Figure  4B.  The  maximum  change  in  Coulomb  stress  is 
taken  from  the  point  on  the  ITSC  where  Coulomb  stress  reaches  a  maximum  value.  The 
results  shown  are  for  calculations  assuming  a  tapered  slip  distribution  along  the  source 
earthquake. 
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Figure  7.  Coulomb  stress  models  for  a  teleseismically  recorded  earthquake  (Mw  =  5.7, 
26  April  2001)  on  the  Siqueiros  transform  system.  Earthquake  location  is  shown  by 
white  star  on  each  panel.  A:  Location  map  shows  the  segmentation  geometry  of 
Siqueiros.  Outlined  region  indicates  area  investigated  in  B,  C,  and  D.  B:  Bathymetric 
map  overlain  by  the  geologic  interpretations  (thin  white  lines)  of  Fomari  et  al.  (1989). 
White  circles  indicate  the  locations  of  the  170  aftershocks.  C:  Calculated  Coulomb  static 
stress  changes  on  secondary  strike-slip  receiver  faults  with  the  same  dip,  strike,  and  rake 
as  the  source  earthquake.  Source  fault  parameters  from  Harvard  CMT  focal  mechanism: 
strike  =  263°,  dip  =  81°,  and  scalar  moment  =  4.49  x  1024  dyne-cm.  We  assume  rupture 
length  =  5  km,  width  =  5  km  and  slip  =  0.36  m.  We  used  a  tapered  slip  distribution  and 
stresses  were  sampled  on  a  horizontal  plane  at  a  depth  of  2  km.  Fox  et  al.  (2001) 
estimates  a  lower  threshold  for  earthquakes  recorded  within  the  hydrophone  array  of  Mb 
=  1.0- 1.8.  The  first  45  aftershocks  (shown  as  white  circles)  occurred  within  10  hours  of 
the  mainshock  and  fall  along  fault  segment  S3.  The  0.15  bar  and  0.05  bar  contours  are 
shown  as  solid  black  lines.  We  observe  that  -31%  of  the  first  45  aftershocks  occurred  in 
regions  with  stress  increases  >0.15  bars,  and  -56%  in  areas  with  stress  increases  >  0.05 
bars.  D:  Calculated  Coulomb  static  stress  changes  on  secondary  normal  faults  dipping 
60°  and  parallel  to  the  ITSC  SB.  The  later  125  aftershocks  (shown  as  white  circles) 
correspond  well  with  predictions  of  increased  normal  stress  thus  suggesting  that  they 
might  be  associated  with  normal  receiver  faults.  We  observe  that  -63%  of  the  latter  125 
aftershocks  occurred  in  regions  with  stress  increases  >0.15  bars,  and  -  90%  occurred  in 
areas  with  stress  increases  >  0.05  bars. 
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Examples  of  Possible  Stress  Interaction 

Hydroacoustic  monitoring  of  the  EPR  (Fox  et  al.,  2001)  has  allowed  us  to  investigate 
several  moderate  sized  earthquakes  to  determine  the  role  of  transform  segmentation  in 
earthquake  processes.  Here,  we  have  chosen  one  Mw  =  5.7  teleseismically  recorded 
earthquake  that  occurred  in  April  2001  along  the  S3  segment  of  the  Siqueiros  transform 
fault,  which  appears  to  have  triggered  seismicity  on  the  S2  segment  as  well  as  on  the 
ITSC  SB  (Fig.  7).  Coulomb  stress  models  were  calculated  using  the  source  earthquake 
focal  mechanism  recorded  in  the  Harvard  CMT  catalog,  and  the  earthquake  source 
location  was  taken  from  the  hydroacoustic  earthquake  catalog. 

The  mainshock  ruptured  the  24-km  long  S3  fault  segment  ~10  km  from  its 
intersection  with  the  8  km-long  ITSC  SB.  Approximately  170  aftershocks  were  recorded 
hydroacoustically  within  24  hr  and  50  km  of  the  epicenter  of  the  mainshock.  The  first  45 
aftershocks  (Fig.  1C)  occurred  along  the  fault  segment  S3  within  10  hr  of  the  mainshock, 
and  were  spatially  truncated  by  the  ITSCs  SB  and  SC.  These  aftershock  locations 
correspond  well  with  areas  of  predicted  static  stress  increase  along  the  source  fault  and  on 
secondary  strike-slip  receiver  faults  along  segment  S3.  Furthermore,  the  termination  of 
the  aftershocks  east  of  the  source  fault  at  ITSC  SC  agrees  with  predicted  estimates  for 
segment  interaction. 

The  subsequent  125  aftershocks  (Fig.  7D)  occurred  on  ITSC  SB  and  the  eastern  end 
of  the  adjacent  fault  segment  S2.  The  static  stress  models  assuming  transform-parallel 
strike-slip  receiver  faults  do  not  predict  this  pattern  of  seismicity  (Fig.  7C).  Geologic 
interpretations  of  Siqueiros  transform  by  Fomari  et  al.  (1989)  indicate  several  ITSC- 
parallel  faults  flanking  ITSC  SB  to  the  east,  west,  and  north  (Fig.  7B).  As  the  latter  125 
aftershock  locations  correspond  well  with  areas  of  increased  Coulomb  stress  for  normal 
receiver  faults,  we  hypothesize  that  these  aftershocks  might  be  associated  with  triggered 
seismicity  on  the  normal  faults  mapped  by  Fornari  et  al.  (1989).  The  ITSC  parallel 
tectonic  fabric  was  created  by  ITSCs  SA  and  SB,  which  appear  to  have  slow-spreading 
ridge  morphology  (Fomari  et  al.,  1989);  this  may  account  for  the  development  of 
seismically  active  normal  faults.  Another  possibility  is  that  these  aftershocks  may  reflect 


38 


a  diking  event  near  ITSC  SB  resulting  from  mainshock-triggered  decreases  in  confining 
pressure.  Dynamic  stress  changes  might  also  trigger  aftershocks,  but  are  difficult  to 
evaluate  due  to  the  lack  of  detailed  mainshock  rupture  models. 

Conclusions 

Detailed  analysis  of  earthquakes  on  transform  systems  at  the  EPR  suggests  that 
segmentation  geometry  plays  an  important  role  in  how  slip  is  accommodated  at  fast¬ 
slipping  oceanic  transforms.  Results  of  Coulomb  stress  models  suggest  that  the  length  of 
the  ITSC  that  offsets  two  transform  fault  strands  will  determine  whether  the  adjacent 
fault  segments  will  interact  by  static  stress  transfer.  If  the  ITSC  is  sufficiently  long,  the 
adjacent  segments  will  be  decoupled  and  behave  independently  of  each  other.  This  is 
particularly  important  in  studies  of  earthquakes  at  oceanic  transforms  since  a  long 
oceanic  transform  system,  could  be  composed  of  several  decoupled  fault  segments. 
Moreover,  we  illustrate  that  the  thermal  effect  of  ITSCs  may  reduce  the  thickness  of  the 
brittle  layer,  thus  decreasing  the  potential  rupture  area  and  the  maximum  seismic  moment 
of  an  oceanic  transform  fault  system.  Finally,  we  suggest  that  transform  earthquakes  may 
have  the  potential  to  trigger  seismicity  on  secondary  normal  faults  flanking  ITSCs. 

Acknowledgements 

This  work  was  supported  by  a  NSF  Graduate  Research  Fellowship  (P.  M.  Gregg), 
NSF  Grant  OCE-0221386  (D.  K.  Smith  and  J.  Lin),  and  the  WHOI  Deep  Ocean 
Exploration  Institute  (J.  Lin).  We  are  grateful  for  helpful  discussions  with  C.  Williams, 
M.  Behn,  M.  Boettcher,  J.  Cann,  D.  Fomari,  D.  Forsyth,  J.  McGuire,  H.  Schouten,  and 
the  WHOI  geophysics  group,  and  to  D.  Bohnenstiehl,  E.A.  Okal,  and  S.  Toda  for 
constructive  and  insightful  reviews. 


39 


References  Cited 


Abercrombie,  R.E.,  and  Ekstrom,  G.,  2001,  Earthquake  slip  on  oceanic  transform  faults: 
Nature,  v.  410,  p.  74-77. 

Boettcher,  M.S.,  and  Jordan,  T.H.,  2004,  Earthquake  scaling  relations  for  mid-ocean 
ridge  transform  faults:  Journal  of  Geophysical  Research-Solid  Earth,  v.  109,  B12, 
doi :  1 0. 1 029/2004J  B003 1 1 0. 

Bohnenstiehl,  D.R.,  Tolstoy,  M.,  and  Chapp,  E.,  2004,  Breaking  into  the  plate:  A  7.6  Mw 
fracture-zone  earthquake  adjacent  to  the  Central  Indian  Ridge:  Geophysical  Research 
Letters,  v.  31,  B2,  doi:  10.1029/2003GL0 18981. 

Bohnenstiehl,  D.R.,  Tolstoy,  M.,  Dziak,  R.P.,  Fox,  C.G.,  and  Smith,  D.K.,  2002, 
Aftershock  sequences  in  the  mid-ocean  ridge  environment:  an  analysis  using 
hydroacoustic  data:  Tectonophysics,  v.  354,  p.  49-70. 

Canales,  J.P.,  Detrick,  R.,  Toomey,  D.R.,  and  Wilcock,  S.D.,  2003,  Segment-scale 
variations  in  the  crustal  structure  of  150-300  kyr  old  fast  spreading  oceanic  crust 
(East  Pacific  Rise,  8°15'N-10°5'N)  from  wide-angle  seismic  refraction  profiles: 
Geophysical  Journal  International,  v.  152,  p.  766-794. 

Dziak,  R.P.,  Fox,  C.G.,  and  Embley,  R.W.,  1991,  Relationship  between  seismicity  and 
geologic  structure  of  the  Blanco  Transform  Zone:  Marine  Geophysical  Researches,  v. 
13,  p.  203-208. 

Fomari,  D.J.,  Gallo,  D.G.,  Edwards,  M.H.,  Madsen,  J.A.,  Perfit,  M.R.,  and  Shor,  A.N., 
1989,  Structure  and  topography  of  the  Siqueiros  transform-fault  system  -  Evidence 
for  the  development  of  intra-transform  spreading  centers:  Marine  Geophysical 
Researches,  v.  1 1,  p.  263-299. 

Forsyth,  D.W.,  Yang,  Y.J.,  Mangriotis,  M.D.,  and  Shen,  Y.,  2003,  Coupled  seismic  slip 
on  adjacent  oceanic  transform  faults:  Geophysical  Research  Letters,  v.  30, 
doi:  10.1029/2002GL01 6454. 


40 


Fox,  C.G.,  Matsumoto,  H.,  and  Lau,  T.-K.A.,  2001,  Monitoring  Pacific  Ocean  Seismicity 
from  an  autonomous  hydrophone  array:  Journal  of  Geophysical  Research,  v.  106,  p. 
4183-4206. 

Goff,  J.A.,  Fomari,  D.J.,  Cochran,  J.R.,  Keeley,  C.,  and  Malinvemo,  A.,  1993,  Wilkes 
transform  system  and  nannoplate:  Geology,  v.  21,  p.  623-626. 

Harris,  R.A.,  and  Day,  S.M.,  1993,  Dynamics  of  fault  interaction  -  Parallel  strike-slip 
faults:  Journal  of  Geophysical  Research-Solid  Earth,  v.  98,  p.  4461-4472. 

Hekinian,  R.,  Bideau,  D.,  Cannat,  M.,  Francheteau,  J.,  and  Hebert,  R.,  1992,  Volcanic 
activity  and  crust  mantle  exposure  in  the  ultrafast  Garrett  transform-fault  near  13- 
degrees-28s  in  the  Pacific:  Earth  and  Planetary  Science  Letters,  v.  108,  p.  259-275. 

King,  G.C.P.,  Stein,  R.S.,  and  Lin,  J.,  1994,  Static  stress  changes  and  the  triggering  of 
earthquakes:  Bulletin  of  the  Seismological  Society  of  America,  v.  84,  p.  935-953. 

Lonsdale,  P.,  1989,  Segmentation  of  the  Pacific-Nazca  Spreading  Center,  1°N-20°S: 
Journal  of  Geophysical  Research,  v.  94,  p.  12,197-12,225. 

McGuire,  J.J.,  Jordan,  T.H.,  and  Lin,  J.,  2002,  Complexities  of  transform  boundaries  in 
the  oceans,  in  Stein,  S.,  and  Freymueller,  J.,  eds.,  Plate  Boundary  Zones,  Volume  30: 
Washington,  D.C.,  AGU,  p.  425. 

McKenzie,  D.P.,  1969,  Speculations  on  the  consequences  and  causes  of  plate  motions: 
Geophysical  Journal  of  the  Royal  Astronomical  Society,  v.  18,  p.  1-32. 

Menard,  H.W.,  and  Atwater,  T.,  1969,  Origin  of  fracture  zone  topography:  Nature 
(London),  v.  222,  p.  1037-1040. 

Perfit,  M.R.,  Fornari,  D.J.,  Ridley,  W.L,  Kirk,  P.D.,  Casey,  J.,  Kastens,  K.A.,  Reynolds, 
J.R.,  Edwards,  M.,  Desonie,  D.,  Shuster,  R.,  and  Paradis,  S.,  1996,  Recent  volcanism 
in  the  Siqueiros  transform  fault:  Picritic  basalts  and  implications  for  MORB  magma 
genesis:  Earth  and  Planetary  Science  Letters,  v.  141,  p.  91-108. 

Pockalny,  R.A.,  1997,  Evidence  of  transpression  along  the  Clipperton  Transform: 
Implications  for  processes  of  plate  boundary  reorganization:  Earth  and  Planetary 


41 


Science  Letters,  v.  146,  p.  449-464. 

Pockalny,  R.A.,  Fox,  P.J.,  Fornari,  D.J.,  Macdonald,  K.C.,  and  Perfit,  M.R.,  1997, 
Tectonic  reconstruction  of  the  Clipperton  and  Siqueiros  Fracture  Zones:  Evidence  and 
consequences  of  plate  motion  change  for  the  last  3  Myr:  Journal  of  Geophysical 
Research-Solid  Earth,  v.  102,  p.  3167-3181. 

Searle,  R.C.,  1983,  Multiple,  closely  spaced  transform  faults  in  fast-slipping  fracture 
zones:  Geology,  v.  1 1,  p.  607-610. 

Stein,  S.,  and  Wysession,  M.,  2003,  An  introduction  to  seismology,  earthquakes,  and 
earth  structure:  Malden,  MA,  Blackwell  Publishing  Ltd.,  498  p. 

Toda,  S.,  and  Stein,  R.S.,  2000,  Did  stress  triggering  cause  the  large  off-fault  aftershocks 
of  the  25  March  1998  M-w=8.1  Antarctic  plate  earthquake?:  Geophysical  Research 
Letters,  v.  27,  p.  2301-2304. 

Toda,  S.,  Stein,  R.S.,  Reasenberg,  P.A.,  Dieterich,  J.H.,  and  Yoshida,  A.,  1998,  Stress 
transferred  by  the  1995  M-w  =  6.9  Kobe,  Japan,  shock:  Effect  on  aftershocks  and 
future  earthquake  probabilities:  Journal  of  Geophysical  Research-Solid  Earth,  v.  103, 
p.  24543-24565. 


42 


Chapter  2:  Spreading  rate  dependence  of  gravity  anomalies  along  oceanic 
transform  faults 


Abstract 

Mid-ocean  ridge  morphology  and  crustal  accretion  are  known  to  depend  on  the 
spreading  rate  of  the  ridge.  Slow-spreading  mid-ocean  ridge  segments  exhibit  significant 
crustal  thinning  towards  transform  and  non-transform  offsets'  '\  which  is  thought  to 
arise  from  a  three-dimensional  process  of  buoyant  mantle  upwelling  and  melt  migration 
focused  beneath  the  centres  of  ridge  segments' 2  4  7‘9'10'12.  In  contrast,  fast-spreading  mid¬ 
ocean  ridges  are  characterized  by  smaller,  segment-scale  variations  in  crustal  thickness, 
which  reflect  more  uniform  mantle  upwelling  beneath  the  ridge  axis' '  L\  Here  we  present 
a  systematic  study  of  the  residual  mantle  Bouguer  gravity  anomaly  of  19  oceanic 
transform  faults  that  reveals  a  strong  correlation  between  gravity  signature  and  spreading 
rate.  Previous  studies  have  shown  that  slow-slipping  transform  faults  are  marked  by  more 
positive  gravity  anomalies  than  their  adjacent  ridge  segments'  2'4  <’,  but  our  analysis 
reveals  that  intermediate  and  fast-slipping  transform  faults  exhibit  more  negative  gravity 
anomalies  than  their  adjacent  ridge  segments.  This  finding  indicates  that  there  is  a  mass 
deficit  at  intermediate  and  fast-slipping  transform  faults,  which  could  reflect  increased 
rock  porosity,  serpentinization  of  mantle  peridotite,  and/or  crustal  thickening.  The  most 
negative  anomalies  correspond  to  topographic  highs  flanking  the  transform  faults,  rather 
than  to  transform  troughs  (where  deformation  is  probably  focused  and  porosity  and 
alteration  are  expected  to  be  greatest),  indicating  that  crustal  thickening  could  be  an 
important  contributor  to  the  negative  gravity  anomalies  observed.  This  finding  in  turn 
suggests  that  three-dimensional  magma  accretion  may  occur  near  intermediate  and  fast¬ 
slipping  transform  faults. 
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LETTERS 


Spreading  rate  dependence  of  gravity  anomalies 
along  oceanic  transform  faults 

Patricia  M.  Gregg1,  Jian  Lin2,  Mark  D.  Behn2  &  Laurent  G.  J.  Montesi2 


Mid- ocean  ridge  morphology  and  crustal  accretion  are  known  to 
depend  on  the  spreading  rate  of  the  ridge.  Slow-spreading  inid- 
ocean- ridge  segments  exhibit  significant  crustal  thinning  tow  ards 
transform  and  non-transform  offsets1'12*  which  is  thought  to  arise 
from  a  three-dimensional  process  of  buoyant  mantle  upwelling 
and  melt  migration  tocuscd  beneath  the  centres  of  ridge  seg- 
inents,'2’4'7-*'1®-‘I.  In  contrast*  fast- spreading  mid-ocean  ridges  are 
characterized  by  smaller,  segment-scale  variations  in  crustal  thick¬ 
ness,  which  reflect  more  uniform  mantle  upwelling  beneath  the 
ridge  axis,vu.  Here  we  present  a  systematic  study  of  the  residual 
mantle  Bouguer  gravity  anomaly  of  19  oceanic  transform  f  aults 
that  reveals  a  strong  correlation  between  gravity  signature  and 
spreading  rate.  Previous  studies  have  shown  that  slow-slipping 
transform  faults  are  marked  by  more  positive  gravity  anomalies 
than  their  adjacent  ridge  segments1*2’4,6*  but  our  analysis  reveals 
that  intermediate  and  last -slipping  transform  faults  exhibit  more 
negative  gravity  anomalies  than  their  adjacent  ridge  segments. 
This  finding  indicates  that  there  is  a  mass  deficit  at  intermcdiale- 
and  fast-slipping  transform  faults,  which  could  reflect  increased 
rock  porosity*  serpent  ini/at  ion  of  mantle  peridotite*  and/or  crus¬ 
tal  thickening.  The  most  negative  anomalies  correspond  to  topo 
graphic  highs  flanking  the  transform  faults*  rather  than  to 
transform  troughs  (where  deformation  is  probably  focused  and 
porosity  and  alteration  are  expected  to  be  greatest  I*  indicating  that 
crustal  thickening  could  he  an  important  contributor  to  the  nega¬ 
tive  gravity  anomalies  observed.  This  finding  in  turn  suggests  that 
three-dimensional  magma  accretion  may  occur  near  inlermediate- 
and  fast-slipping  transform  faults. 

We  analyse  the  residual  mantle  Bouguer  gravity  anomaly  (R MBA) 
from  19  oceanic  transform  fault  systems  and  their  adjacent  ridge 
segments  at  ultraslaw  to  fast  spreading  mid -ocean  ridges  (Supple¬ 
mentary  Fig.  1)  with  fault  slip  rates  of  13-14.5cmyr~ 1  and  fault 
lengths  of  70-350 km  (see  Methods).  RMBA  is  used  to  infer  local 
density  variations  relative  to  a  reference  model  that  inciudcsa  6-km 
thick  crust.  A  negative  RMBA  indicates  a  mass  deficit,  which  can  be 
due  to  thickened  crust,  increased  porosity,  and/or  serpentini/ed 
mantle.  Conversely,  a  positive  RMBA  typically  implies  crustal  thin¬ 
ning.  Our  results  reveal  a  systematic  disparity  in  the  gravity  signature 
between  ultra  slow  /slow  slipping  transform  faults  and  intermediate  / 
fast  slipping  transform  faults. 

To  illustrate  these  differences,  we  compare  the  RMBA  at  the  fast 
slipping  Siqueiros  transform  to  the  RMBA  at  the  slow -slipping 
Atlantis  transform  (Fig.  I).  The  I  SO- km  long  Siqueiros  transform 
located  at  -'-8.3'N  on  the  Bast  Pacific  Rise  is  slipping  at 
*•**■  1 1.8cm  yr  1  and  includes  five  fault  segments  connected  by  four 
intra- transform  spreading  centres14  (Fig.  la).  The  ridge  segments 
to  the  north  and  south  of  Siqueiros  are  typical  of  fast -spreading 
ridges  with  axial  highs  as  shallow  as  2.0-2.8  km.  l*he  ridge  transform 


intersections  bounding  Siqueiros  display  overshooting  ridge  tips16 
typical  of  ridge  transform  intersections  at  magrriatically  robust  ridge 
segments  ( Supplementary  Figs  2-5).  In  contrast,  the  70  km  king 
Atlantis  transform  at  30"  N  on  the  Mkl  Atlantic  Ridge  is  slipping  at 
*~'2.4cmyr~ 1  and  contains  a  single  strike  slip  fault  zone  (Fig  lb). 
The  ridge  segments  adjacent  to  the  Adantis  transform  display  slow- 
spreading  rift  valleys  and  a  wide  range  of  axial  depthv 

As  seen  in  previous  gravity  studies  of  the  Atlantis  transform2*717, 
more  negative  RMBA  is  associated  with  adjacent  ridge  segment  ten 
tres  (up  to  ^80mC»al  less  than  the  transform  offset*,  whereas  the 
most  positive  RMBA  is  kkrafed  in  the  transform  fault  and  fracture 
zone  (Fig.  Id).  In  contrast*  at  Siqueiros  more  positive  RMBA  is 
observed  along  the  ad  jacent  ridge  segments  and  more  negative  values 
are  found  within  the  transform  fault  and  fracture  zone  domains  (up 
to  —40mCiaJ  less  than  the  adjacent  rklge  segments).  RMBA  mini 
mums  often  correspond  to  bathymetric  highs  ak>rtg  the  eastern  ami 
southern  flanks  of  the  fault  (  Fig.  1c). 

The  difference  in  t  he  average  RM  BA  values  between  transform  fault 
domains  and  their  adjacent  ridge  segments — ARMBAT  R — varies 
systematical!)'  with  spreading  rate  for  t  Ik*  1 9  transform  faults  included 
in  this  study  (Fig.  2.  sec  Methods  section  for  ARMBAt  K  cakula 
lion).  ARMBAt  r  >  f)  implies  more  positive  RMBA  akmg  the1  trans¬ 
form  fault  domain  than  along  the  adjacent  ridge  segments,  whereas 
ARMBAt  r  <0  implies  more  negative  values  within  the  transform 
fault  domain.  Along  ultrasknv  and  slow  spreading  ridges,  the  average 
RMBA  in  transform  fault  domains  is  systematically  more  poshivethan 
along  the  adjacent  ridge  segments  with  ARMBAt  r  of  up  to  45  rn<  ial. 
However,  for  spreading  rates  exceeding  — 5cmyr  ‘  the  average 
RMBA  within  transform  fault  domains  becomes  more  negative  than 
on  the  adjacent  ridge  segments.  Most  intermediate  and  fast  slipping 
transforms  are  associated  with  ARMBAt  r  of  —4  to  —  I  Initial, 
whereas  the  Blanco  transform  is  associated  with  a  ARMBAt  r  °f 
—  40m(ial  (Fig.  2  and  Supplementary'  Table  1). 

Previous  studies  attribute  positive  RMBA  along  slow -slipping 
transform  faults  to  decreased  crustal  accretion  towards  the  end  of 
ridge  segments  and  therefore  crustal  thinning  within  the  transform 
fault  domain1-2,4*.  If  the  negative  RMBA  observed  at  intermediate-  and 
fast -slipping  transforms  also  reflects  variations  in  crustal  thickness, 
these  anomalies  indicate  crustal  thickening  within  the  transform 
domain.  However,  negative  RMBA  can  also  arise  from  other  factors, 
such  as  cracking  and  alteration  at  transform  fault  offsets5'*'10.  There 
fore,  we  use  two-dimensional  (2D)  gravity  models  to  estimate  the 
amplitude  of  ARMBAt  r  corresponding  to:  ( I  )  increased  rock  poro 
sity,  (2)  serpent  inizat  ion  of  mantle  peridotite:  and  (3)  crustal  thick 
cning  (  Fig.  3,  Supplementary  Fig.  7,  and  Supplementary  Methods). 

Low  seismic  velocities  recorded  in  the  Qipperton  transform  and 
near  the  Siqueiros  transform  have  been  attributed  to  increased  crack¬ 
ing  and  porosity  within  the  transform  domain151*.  The  observed 
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Figure  1  Comparison  of  bathymetry  and  RMBA  of  the  Siqueiros  transform 
on  thK?  East  Pacific  Rise  and  Atlantis  transform  on  the  Mid-Atlantic  Ridge, 
at  the  same  map  scale,  a.  Bathymetric  map  lor  the  Siqueiros  transform  with 
ideologic  interpretation**.  Solid  black  lines  indicate  Lhe  seafloor! ahric,  circles 
show  locations  of  seamounts,  and  dashed  black  lines  show  the  locations  ol 
the  fracture  /ones  and  trail  storm  tault.  The  solid  white  line  indicates  the 
location  of  the  plate  boundary  used  in  the  3D  thermal  model.  II  SC  Is  are 
labelled  A,  B,  ('  and  D*.  NV/,  neo« volcanic  /one;  OSC,  overlapping 


Figure  2  Compilation  of  ARMBAt_r  values  for  the  19  transform  systems 
analysed.  The  new  RMBA  calculations  from  thus  study  are  marked  with  solid 
circles.  Previous  RMBA  studies  are  indicated  by  the  open  symbols:  squared 
triangle*2.  crtKSJ0,  star*1,  diamond'.  Abbreviations  for  ridges:  SWIR. 
Southwest  1  ndian  Ridge; MAR.  Mid  Atlantic  Ridge;  I  dl-  K.  luan  de  1- ucu  Ridge. 
SE1R,  Southeast  Indian  Ridge;  EPR.  East  Pacific  Rise.  Abbreviations  for 
transform  faults:  Bo.  Bouvet;  IO,  Idas  Orvadas;  A  B,  Andrew  Bain;  PE,  Prince 
Edward;  All,  Atlantis  11;  Mo,  Moshesh;  Ma,  Marion;  At,  Atlantis;  11s,  Hayes 
Kn,  Kane;  FT,  Fifteen  Twenty;  AFZ,  Ascension;  RCi,  Rio  Grande;  Bl,  Blanco; 
S2.  SL1R2;  SI.  SUR1;  Q,  Clippertun;  Sq.  Siqueiros;  Wi,  Wilkes;  Cia.  Ciarrett. 


spreading  centre,  b  Bathymetric  map  tor  the  Atlantis  transform  The  solid 
black  line  indicates  location  ol  the  plate  boundary  used  in  the  3D  thermal 
model,  and  the  dashed  portion  indicates  the  Atlantis  transform  lault 
c,  Calculated  RMBA  map  for  the  Siqueiros  transform  lault  with  geologic 
interpretation  overlaid  from  a.  d,  Calculated  RMBA  map  lor  the  Atlantis 
transform  fault.  Seabeam  bathymetry  compiled  by  the  RIDCiE  multibeam 
synthesis  profe*.!  (hllps/'occan  ndge.ldeo. Columbia. edu'generaL'htmL 
home  Itttnl !. 

decrease  in  seismic  velocities  of  I  km  s  1  in  the  (  lipperton  trans 
form  fault1*  can  be  explained  by  porosity  increases  in  the  range  of 
I  3.5%  for  crack  aspect  ratios  of  0.01  -0.1  « ref*.  19)  (Calculations  in 
Kig.  3  demonstrate  that  an  average  porosity  increase  of  1  -7%  within 
the  upper  crust  down  to  the  brittle  ducule  transition  uOO  (!'•  is 
required  to  produce  the  negative  ARMBAj  k  values  of  up  to 
—  1 -Initial  observed  at  fast  slipping  transforms  <  1-ig.  2'  However, 


I  *”‘1 
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Figure  3  i  Results  of  2D  forward  models  showing  the  predicted 
ARMBAt  _  «•  Data  are  shown  as  a  function  of  increase  in  crustal  thickness 
through  layer  2A  i  dotted  black  linel  or  overall  crustal  thickening  (solid  black 
linel;  increases  in  porosity  10,  red  line);  or  the  presence  of  serpentini/ed 
mantle  i  green  lines  indicate  percentage;  green  shading  indicates  serpentine 
stability  field  .  Ah  indicates  change  in  layer  thickness.  Grey  shading  indicates 
the  observed  ARMBAr  -  r  at  last  slipping  transforms.  The  kink  in  the  curve 
for  increase  in  layer  2A  corresponds  to  thickening  ol  layer  2A  beyond  tin.* 
thickness  of  layer  2B  in  the  reference  model.  See  Supplementary  Fig.  7  lor 
model  details. 
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Figure  4|  Lateral  variations  in  crustal  thickness  required  to  explain  the 
observed  RMBAinthe  Siqueiros  transform  system  (Fig.  1c).  Excess  crust  is 
defined  as  the  deviation  from  the  reference  crustal  thickness  of  tv  km, 
assuming  sea  water,  crust  and  mantle  densities  of  1,030,  2,730  and 
3,300  kg m  5.  respectively. 


an  increase  in  porosity  in  excess  of  5%  at  slow-slipping  transforms 
will  negate  the  effect  of  crustal  thinning  and  result  in  negative 
ARMBAt_r  •  Supplementary  Fig.  7).  Consequently,  a  porosity 
increase  of  2-5%  is  consistent  with  observed  ARM  15 AT  _  R  and  seis¬ 
mic  data  Jt  both  slow  and  fast  slipping  transform  faults. 

Serpentinization  of  mantle  peridotitc  can  also  generate  negative 
ARM  BA-i  -  r.  I  low  ever,  the  calculated  500  C  isotherm,  which  repre¬ 
sents  die  upper  limit  of  serpentine  st al»il it  v20,  is  located  within  die 
crust  at  all  last  slipping  transform  faults  considered  in  this  study 
(see  Methods)  and  therefore  predicts  no  mantle  serpentinization. 
Invoking  hydrothermal  cooling  (see  Supplementary  Methods) 
depresses  the  500  0  isotherm  to  an  average  depth  of  6.4  km  at 
Siqueiros  transform  fault.  Assuming  100%  serpentinization  (density 
of  2,550  kg  m~3)  of  the  mantle  where  temperatures  are  <300  :C 
results  in  a  negative  ARMBAr  _  R  of  0  to  —  lOmCial,  which  is  slighdv 
less  than  the  values  of  up  to  —  1 4  mCial  observed  at  fast  slipping  trans¬ 
form  faults  ( Fig.  3).  However,  if  100%  serpentinization  were  pervasive 
above  the  500  *C  isotherm  at  slow  slipping  transform  faults,  it 
would  negate  the  effects  of  crustal  thinning  and  produce  negative 
ARM  BA,  R<  —  lOOmfial  (Supplementary  Fig  7g).  Therefore, 
although  low  degrees  of  mantle  serpentinization  might  contribute 
to  the  observed  negative  RMBA  along  fast  slipping  transform  faults, 
current  thermal  models  for  oceanic  transform  faults  make  it  unlikely 
that  serpentinization  alone  can  explain  the  systematic  differences 
between  fast-  and  slow-slipping  transforms. 


Finally,  we  estimate  the  magnitude  of  variations  in  crustal  thick 
ness  required  to  explain  the  RMBA  observed  along  the  Siqueiros 
transform  by  inverting  the  observed  RMBA  relative  to  a  reference 
model  with  average  crustal  density  of  2,730 kg  m-3  (ref.  1 )  (Fig.  4). 
Along  the  eastern  domain  of  Siqueiros  transform,  we  predict  an 
excess  crustal  thickness  of  >1.5  km,  while  in  the  western  domain  of 
the  transform  and  in  the  off  axis  fracture  zones  excess  crust  ranges 
from  0.5- 1.5  km.  The  ARMBAj  .  R  observed  at  fast  slipping  trans 
form  faults  can  also  be  produced  by  increasing  laver  2A  bv  0.25- 
0.6  km  (Fig.  3). 

Although  increased  rock  porosity  and  mantle  serpentinization 
could  both  contribute  to  the  observed  negative  RMBA  along  inter 
mediate-  and  fast-slipping  transform  faults,  we  believe  that  local 
crustal  thickening  is  probably  the  dominant  mechanism.  Areas  of 
localized  deformation,  such  as  fault  zones,  are  expected  to  exhibit 
enhanced  cracking  and  alteration.  However,  the  most  negative 
RMBA  values  observed  at  intermediate-  and  fast  slipping  transform 
faults  correspond  to  bathymetric  highs  flanking  the  transforms 
instead  of  the  transform  troughs  where  deformation  is  likely  to  be 
localized.  Local  excess  crust  within  intermediate-  jnd  fast -slipping 
transform  fault  domains  mav  result  from  a  combination  of  magmatic- 
accretion  mechanisms  including:  ( 1 )  lateral  transport  of  excess 
magma  from  the  ridge  axis  into  the  transform  fault  domain16-21;  (2) 
magmatic  accretion  at  intra  transform  spreading  centres1*'23;  (3)  the 
pooling  of  lava  within  topographic  lows2';  and/or  (4i  leaky'  mag¬ 
matic  accretion  along  the  entire  transform  fault24. 

Lateral  dyke  propagation  from  a  central  magma  chamber  has 
been  proposed  to  explain  crustal  thickness  variations  along  slow- 
spreading  ridge  segments25.  At  magmatically  robust  segments  along 
intermediate  and  fast  spreading  ridges,  a  continuous  magma  cham 
her  may  exist  beneath  the  entire  segment  and  feed  dvkes  from  the 
ridge  transform  intersect k>n  into  the  transform  fault  domain.  A  pos¬ 
sible  geologic  expression  of  this  mechanism  is  overshooting  ridge 
tips,  which  are  common  bathymetric  features  of  robust  ridge  seg 
menls  at  intermediate-  and  fast- spreading  ridge  transform  intersec¬ 
tions16  (see  Supplementary  Figs  2-5).  Bathymetry  and  RMBA  of 
overshooting  ridge  tips  at  intermediate  and  fast -spreading  ridges 
suggest  that  dykes  propagating  along  the  ridge  axis  penetrate  past 
the  transform  fault  into  juxtaposed  old  oceanic  crust,  curve  in  the 
direction  of  the  transform,  and  locally  increase  the  crustal  thickness 
of  an  already  full  thickness  crust  (Fig.  5bc  Geologic  evidence  from 
the  northwestern  segment  of  Blanco  transform  indicates  significant 


Figure  5  |  A  spreading-rate-depeodent  model  of  crustal  accretion  and 
mantle  up  welling  based  on  observed  RMBA  calculations  and  morphological 
features  at  transform  fault  systems  on  slow-  and  fast-spreading  ridges. 

Large  yellow  arrows  indicate  plate  motion  vectors;  smaller  black  arrows 
indicate  mantle  flow.  The  grey  area  with  white  outline  represents  the 
lithosphere  as  defined  by  the  700  T  isotherm,  while  ihe  crust  is  identified  by 
dotted  pattern.  The  blue  cracks  and  the  green  wiggles  indicate  regions  of 
increased  porosity  and  serpentinization,  respectively.  *.  Slow-spreading 


ridges  exhibit  focused  mantle  upweiling  at  segment  centres  and  crustal 
thinning  towards  segment  ends  and  transform  faults,  b,  l  ast-spreading 
ridges  have  relatively  uniform  mantle  upwel ling. The  negative  ARMBAT_  R 
at  fast-spreading  ridges  may  result  from:  (I)  excess  magma  supply  on  ridge 
segments  causing  over-shooting  ridges  and  dyke  injection  into  the  transform 
fault  domain;  (2)  increased  porosity  within  the  transform;  (3)  actively 
accreting  ITSCs;  (41  pooling  of  extruded  lavas  within  the  transform  fault. 
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increases  in  the  extrusive  lava  layer  with  a  possible  overshooting  ridge 
tip  source"'  I  lowever,  for  this  model  to  explain  the  RMBA  along  the 
entire  transform  fault  (rather  than  111st  in  the  fracture  zone),  trans 
form  parallel  dyke  propagation  must  extend  jt  least  h alt  the  distance 
between  ridge  transform  intersections  (—50-75  km).  Bathymetric 
data  demonstrate  that  overshooting  ridge  tips  wrap  around  and  into 
transform  faults*  but  the  cause  of  this  curvature  is  unclear.  Possible 
explanations  include  response  to  the  local  stress  field  at  the  ridge 
transform  intersection,  ar.cL'or  dyke  propagation  along  preferential 
pathways  created  by  pre  existing  transform  parallel  cracks. 

Hast  slipping  transform  faults  typically  contain  multiple  fault  seg 
merit*  separated  by  infra  transform  spreading  centres  (ITSO11*26 
Several  ITSlisat  Hast  Pacific  Rise  transform  faults  have  been  sampled 
and  contain  fresh  basaltic  lavas  indicating  crustal  accretion  within 
the  transform  fault  domain16*22.  In  some  cases,  such  as  the  tiarrett 
transform  fault  (Supplementary  Fig.  2),  the  TTSC's  iKturin  pull  apart 
basins  and  do  not  appear  to  increase  the  crustal  thickness.  Alter 
natively,  along  the  eastern  portion  of  Siqueiros  transform  (Fig.  I  a • 
ITSCs  are  associated  with  negative  RMBA.  Thus,  ITSCs  mav  thicken 
crust  in  some  localities. 

A.  third  possible  mechanism  for  local  crustal  thickening  is  the 
pooling  of  erupted  lavas  in  topographic  lows  within  the  transform 
fault  domain21  (Fig.  5b).  However,  the  majority  of  the  negative 
RMBA  values  we  have  observed  correspond  with  topographic 
highs,  indicating  that  this  might  not  be  a  dominant  mechanism 
finally,  because  off-axis  melt  is  considered  abundant  at  fast 
spreading  ridges27,  leaky’  magmatic  accretion  might  occur  along 
the  entire  length  of  some  transform  faults24.  Phis  mechanism  could 
i  aiive  increases  in  crustal  thickness  anil  may  lx*  enhanced  by  trans 
tensional  forces  acting  upon  a  transform  fault  owing  to  plate  motion 
reorganizations24. 

Thus  the  observed  negative  RMBA  vjlues  at  intermediate  and 
fast  slipping  transform  faults  may  be  the  result  of  a  combination  of 
mechanisms,  including  increases  in  rock  porosity,  mantle  serpentm 
i/atior.  and/or  crustal  thickening.  Further  detailed  constraints  jre 
required  to  distinguish  among  the  potential  contributors  to  the 
observed  gravity  anomalies.  However,  the  strong  correlation  of  the 
negative  RMBA  with  topographic  highs  flanking  the  transforms  sug 
gests  crustal  thickening  as  a  potentially  important  mechanism  for  at 
least  some  of  the  observed  negative  RMBA  values.  Systematic  occur 
rence  of  crustal  thickening,  it  confirmed,  would  require  a  modified 
rion  of  our  current  understanding  of  how  three-dimensional  magma 
accretion  occurs  along  intermediate-  and  fast-spreading  ridges  in  the 
vicinity  of  transform  fault  discontinuities. 

METHODS  SUMMARY 

Tlw-  imiillr  Bouguer  anomaly  was  ,  .du. Uteri  following  lit?  method  of  ref.  1. 
This  technique  removes  ihe  effect  of  seatloor  topography  and  hkm  of  reterence 
crust  from  the  salellitc  tree  air  gravity^4  The  RMBA  wav  calculated  by  removing 
trom  the  mantle  Bouguer  anomaly  the  effects  of  lithospheric  cooling  as  estimated 
from  a  three  dimensional  (3D)  passive  mantle  upwelling  minlel  calculated 
using  iltr  commercial  finite-element  modelling  software  pat.  luge  CflMSOI. 
M ULTIPHYSICS  3.2  (for  model  details  see  Supplementary  Tables  1  and  2). 
The  RMBA  retleels  the  deviation  of  sub- seafloor  density  structure  trom  the 
retererxe  model  described  in  Supplementary  Table  2.  When  constructing  the 
thermal  models  for  segmented  transform  faults, we  incorporated  thebest  known 
geometry  of  ti  ansform  fault  segments  and  assumed  steady  stale  crustal  accretion 
at  the  ITSCs. 

We  detine  ARMBA,  _  R  *  ((RMBA, -RMBAH, !  *  ( RMBAf- RMBAm )| 72, 
where  RMBAW  and  RM  BA^  are  the  average  RMBA  valueson  the  adjacent  ridge 
segments  1  and  2,  respectively,  and  RMBAj  is  the  average  RM  BA  value  along  the 
transform  fault  domain  (see  Supplementary  Fig.  8).  Where  digital  data  are  avail¬ 
able14**,  RMBAK1  and  RMBAR2  were  calculated  by  averaging  the  values  at  ihe 
ccnlreoleach  adjacenl  ridge  segment  withina  20-km-widcand  50-km-longbox. 
At  the  East  Pacific  Rise,  luan  Je  Fuca  Ridge,  and  Southeast  Indian  Ridge,  where 
the  translorm  tault  zones  are  wider  owing  to  segmentation,  RMBA(  was  calcu¬ 
lated  by  averaging  over  a  20- km- wide  swath  centred  20  km  fnun  the  ridge  trans¬ 
lorm  intersections.  At  the  Mid-Allantic  Ridge  and  Southwest  Indian  Ridge  the 
RMBA  values  were  averaged  over  a  10-km-widr  swath  centred  20  km  fn>m  the 


ridge  transform  intersections.  Where  digital  data  are  unavailable*'  *\  average 
values  were  estimated  trom  published  proliles  and  maps  ot  the  translorm  tadts 
and  (lieir  adjacent  ridge  segments. 

Full  Methods  and  any  associated  references  are  available  in  Ihe  online  ve-sion  J 
Ihe  pacer  at  www.nature.cjm/rialure 
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METHODS 

Because  there  is  an  in  Keren!  non-uniqueness  in  gravity  calculations,  we  have 
..realed  several  2D  tor  ward  miHlels11  to  quantity  Lbe  effects  of  various  median 
isms  on  the  gravity  signature  of  a  ridge-transform  system  (Supplementary  Fig. 
7}.  The  light  blue  arxl  orange  shaded  regions  in  Supplementary  Fig.  7  indicate 
the  extent  ot  ARMBA,  _K  observed  along  transform  faults  at  slow  and  tasi- 
spreaiiing  ridges,  respectively.  The  density  used  tor  each  ot  the  layers  is  indicated 
in  Supplementary  Fig.  7.  These  values  do  vary  in  nature'1-1*. 

In  the  tirst  suite  ot  models  we  varied  the  crustal  thickness  by  A/r  either  entirely 
as  a  change  ot  layer  2A.  layer  2B,  or  proportionally  throughout  the  entire  crust 

•  Supplementary  Figs  7a,  b  and  c.  respectively:.  All  three  mechanisms  act 
to  decrease  ARMBA,  _K  within  the  range  of  fast  slipping  transforms  with 

•  -  1  km  otexcess  material.  Increasing  layer  2Arequircstbe  lexsl  amount  ot  excess 
material  with  0.25-0.6  km.  The  kink  in  model  A  results  trom  increasing  layer  2A 
beyond  the  thickness  ot  layer  2B. 

The  second  suite  ol  mi  dels  considers  increased  ponssity  in  llie  brittle  layer 

•  temperatures  60U  Cl.  Two  3D  thermal  models  were  calculated  to  obtain  the 

average  depth  u>  the  h(J0  C  isotherm  along  the  Siqueiros  transform  fault:  ( :  :  l he 
initial  3D  passive  mantle  upwelling  model,  outlined  in  the  Methods!  the  red  dashed 
line  at  —3km  depth  on  Supplementary  Figs  7d  and  e  and  (2)  a  hydrolhermallv 
ci>olcd  model  in  which  licat  conductivity*  Ls  increased  tenfold  to  4  km  depth  (live 
blue  dashed  line  at  b.4  km  depth  on  Supplementary  Figs7dand  e:.  Density  changes 
due  to  variations  in  porosity  are  defined  as:  pi  «01  =  (1  where  '/> 

is  the  porosity,  ls  the  density  ot  the  host  rock,  arxl  is  the  density  ol  sea 
water.  Assuming  that  the  bril tie-ductile  transition  is  at  -3  km,  porosity  would 
need  to  Ivr  increased  2-7%  to  achieve  the  ARMBA  t  _  R  values  of  magnitude 

14  mGal  observed  at  last  slipping  transform  faults. 

Similar  tabulations  were  performed  h>r  slow-slipping  transforms,  assuming 
a  thicker  brittle  layer,  as  was  calculated  for  the  Atlantis  transform  Hull.  These 
results  iixlieate  that  it  is  difficult  to  exj’lain  the  ARMRA,  R  at  both  fast  and 
slow-slipping  transtorms  simultaneously,  assuming  similar  porosity  structures 
in  both  regions.  Becauseil  isunlikelv  that  there  will  lx:  large  increases  in  porosity 
down  to  Y  km  depth  at  slow- slipping  l  ran  storm  taulls,  even  though  it  is  within 
Lhe  pressure  constraints  suggested  by  ret.  36,  we  have  rerun  the  2D  slow  slipping 
i  ran  storm  model  with  increased  |mrosity  to  fi  km  depth  While  port  wily  i  ix  rea  ses 

in  excess  of  8%  still  act  to  negate  the  effect  of  crustal  thinning  at  slow  slipping 
transform  taulls,  moderate  increases  of  2-5%  are  consistent  with  ARMBA,  _K 
values  at  both  slow-  and  tast-slipping  transform  taulls  (Supplementary  I'tg.  7e). 

The  tinaJ  set  of  forward  models  includes  serpenlinization  iff  the  mantle 
Ixrneath  I  he  transform  fault  (^up|demenl.ary  Figs  TfamlgJ.  As  in  the  porosity 
models,  we  have  incorporated  two  3D  thermal  calculations  at  the  Siqueiros 
transform  to  obtain  average  depths  to  the  SI XI  f.  isotherm,  which  limits  the 
serpentine  stability  lield*'0  to  2. 5km  for  uncooled  and  n.4km  tor  hydrolhermallv 
cixded  lithosphere.  Without  hydrothermal  tooling  lhe50fl  C  Isotherm  lies  well 
within  the  upper  crust  and  there  is  no  predicted  mantle  serpen tiniaa lion. 
However,  if  hydrothermal  axrling  Ls  applied  from  0  to  4  km  depth,  an  end 
member  model  of  100%  serpentin nation  will  result  in  a  ARMBAt-k  of 

lOrnCial,  which  ls  the  average  ot  the  observed  ARMBA,  _  K  lor  tast  slipping 
transform  taulls.  It  a  similar  model  for  mantle  serpent  in  i/at  ion  Ls  applied  at  a 
slow  slipping  transform  fault,  mantle  serpent  in  i/at  ion  down  to  the  hydrolher 
mails  cooled  500  ‘C  isotherm  would  prixluce  such  a  large  negative  RMBA  that  it 
would  completely  offset  the  effect  of  crustal  thinning,  which  is  inconsistent  with 
observations. 

Additionally,  we  calculated  a  modified  RMBA  at  the  Siqueiros  transform 
using  the  hydrolhermallv  doled  thermal  model  and  the  density  for  100%  ser 
pen  tine  (2,550  kg m-’)  in  all  regions  of  the  mantle  I  >6  km  depth;  where  the 
temperature  U  <500  C.  Even  though  this  correction  maximizes  the  effect  of 
serpentine,  it  dec  reaves  the  average  ARMBA  i  _  K  within  the  Siqueiros  transform 
fault  domain  by  only  <2 mGal,  and  reduces  the  calculated  excess  crustal  thick 
ness  by  <200m. 

laboratory  experiments  suggest  that,  if  the  crust  is  thinnedarxl  mantle  upwel¬ 
ling  occurs,  hydrothermal  tluid  penetration  through  thermal  cracking  arxl  grain 
boundary  welting  will  reach  only  —4-6  km  depth*7.  Therefore,  it  is  unlikely  that 
hydrothermal  fluids  wdl  penetrate  to  great  enough  depths  to  serpenlinize  large 
regions  of  mantle  beneath  last-slipping  transform  faults.  Moreover,  depressing 
the  500  C  isotherm  through  hydrothertiul  cooling  is  inconsistent  with  calcu¬ 
lated  earthquake  depths  at  oceanic  transform  faults*4.  On  the  other  hand,  slow 
slipping  transform  faults  display  significant  crustal  thinning  and  are  cooler 
than  intermediate  arxl  fast  slipping  transform  faults,  implying  that  a  larger 
portion  i>f  the  iirxlcrlying  mantle  resides  in  the  serpentine  stability  field. 
However,  the  locally  positive  RMBA  observed  at  slow  slipping  transform  faults 
reflects  crustal  thinning  despite  tlx*  potential  for  significant  serpen  tini/rd  peri  - 
dotite'.  fogethcr.  these  results  indicate  that,  whde  mantle  serpen  tin  ization  may 


ix.ur  al  intermediate  arxl  tast  slipping  l  ran  storm  taulls,  (lie  direct  contribution 

ol  serpent  ini /ali  or.  to  lhe  gravity  signature  Ls  likely  In  lx:  negligible. 
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SI.  Supplementary  Figures 
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Supplementary  Figure  1.  Bathymetric  map  indicating  the  location  of  the  19  transform  faults 
analyzed.  Thin  black  lines  show  the  location  of  plate  boundaries  for  the  areas  of  interest.  New 
residual  mantle  Bougucr  anomaly  (RMBA)  calculations  have  been  performed  for  transform  faults 
in  orange  boxes  (presented  in  Supplemental  Figs.  2-5).  Regions  where  RMBA  calculations  have 
been  previously  published10  are  indicated  by  transform  faults  in  white  boxes.  We  included 
previous  RMBA  studies  only  if  both  transform  fault  and  corresponding  ridges  were  analyzed  by 
the  same  authors  using  the  same  parameters.  All  other  RMBA  studies  we  considered  from  slow- 
spreading  ridges  implied  ARMBAt.r  >  0.  MAR  -  Mid  Atlantic  Ridge,  EPR  -  East  Pacific  Rise, 
JdFR  -  Juan  de  Fuca  Ridge,  SWIR  -  Southwest  Indian  Ridge,  SE1R  -  Southeast  Indian  Ridge, 
AB  -  Andrew  Bain  FZ,  Bo  -  Bouvet  FZ,  Mo  -  Moshcsh  FZ,  PE  -  Prince  Edward  FZ,  All  - 
Atlantis  II  FZ,  10  -  Islas  Orcadas  FZ,  Ma-  Marion  FZ,  SI  -  SEIR1,  S2  -  SEIR2. 
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Supplementary  Figure  2.  (a)  Seafloor  depth  and  (b)  RMBA  for  the  Wilkes  and  Garrett 
transform  faults  on  the  East  Pacific  Rise.  The  solid  black  line  indicates  the  plate  boundary  used 
for  3D  thermal  models  at  each  region6'  .  Arrows  specify  the  location  of  overshooting  ridge  tips 
(ORTs)  at  ridge  transform  intersections  (RTI).  The  doited  line  outlines  the  Wilkes  nanoplate' 
All  supplemental  maps  (Supplemental  Figs  2-6)  arc  at  the  same  map  scale. 


51 


1 1  N 


105’W  103"W  10TW 


Supplementary  Figure  3.  (a)  Seafloor  depth  and  (b)  RMBA  for  the  Siqueiros  and  Clipperton 
transform  faults  on  the  EPR.  The  solid  black  line  indicates  the  plate  boundary  used  for  the  3D 
thermal  models,  and  black  arrows  indicate  location  of  ORTs's.  Note  the  Clipperton  seamount 
chain  to  the  south  and  west  of  the  western  RTI. 
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Supplementary  Figure  4.  (a)  Seafloor  depth  and  (b)  RMBA  for  Southeast  Indian  Ridge 
transform  1  and  transform  2.  The  solid  black  line  indicates  the  plate  boundary  used  for  3D 
thermal  models  at  each  region  ’ H>,  and  black  arrows  indicate  ORTs. 
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Supplementary  Figure  5.  (a)  Seafloor  depth  and  (b)  RMBA  for  the  Blanco  transform  fault  on 
the  Juan  de  Fuca  Ridge.  The  solid  black  line  indicates  the  plate  boundary  used  for  3D  thermal 
models  at  each  region".  The  Blanco  transform  fault  is  the  most  anomalous  of  the  ARMBAt.r 
calculated  (main  text  Fig.  2).  There  are  several  complications  in  this  region,  which  are  not 
incorporated  into  the  3D  thermal  model,  including  the  close  proximity  to  the  Cascadia  subduction 
zone  as  well  as  a  proposed  serpentine  intrusion  forming  Blanco  Ridge  along  the  southeastern 
portion  of  the  transform  fault  (~128°W,  43.5°N)1U2.  While  a  serpentine  block  has  been  proposed 
along  the  southeastern  portion  of  the  transform,  the  most  negative  RMBA  values  are  associated 
with  overshooting  ridge  tip  and  bathymetric  highs  along  the  northwestern  portion  of  the  transform 
fault  (~T30°W,  44.25°N),  and  are  likely  due  to  crustal  increases.  Karson  et  al.  (2002)1’  note  that 
the  extrusive  volcanic  layer  along  the  northwestern  portion  of  Blanco  is  substantially  thicker  than 
what  is  predicted  by  seismic  studies  in  this  region.  An  increase  in  layer  2A  agrees  with  our 
observation  of  negative  RMBA  along  this  portion  of  the  Blanco  transform. 
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Supplementary  Figure  6.  (a)  Seafloor  depth  and  (b)  RMBA  for  the  Hayes  and  Kane  transform 
faults  on  the  Mid-Atlantic  Ridge.  Clear  RMBA  “bulls-cye”  minima  arc  observed  at  segment 
centers,  in  agreement  with  previous  studies  at  slow-spreading  mid-ocean  ridges. 
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Supplementary  Figure  7.  Three  suites  of  forward  models.  The  calculated  ARMBAtr  arc 
plotted  on  the  top  panels.  The  light  blue  and  orange  shaded  region  indicates  the  extent  of 
ARMBAi  r  observed  at  slow  and  fast-slipping  transform  faults,  respectively.  A)  Change  in  layer 
2A  thickness  (dotted  line),  B)  Change  in  layer  2B  thickness  (dashed  line),  C)  Proportional  change 
in  crustal  thickness  (solid  black  line).  D)  Porosity  increase  in  the  brittle  layer  at  an  EPR-type 
fast-spreading  ridge-transform  system.  The  model  illustration  demonstrates  that  porosity  is 
increased  down  to  the  brittle-ductile  transition,  the  600°C  isotherm,  for  both  an  un-coolcd  (red 
dashed  line)  and  hydrothcrmally  cooled  (blue  dashed  line)  thermal  model.  The  red  line  with  solid 
red  triangles  indicates  un-cooled  600°C  isotherm  and  blue  line  with  solid  blue  triangles  indicates 
hydrothermal ly  cooled  600°C  isotherm  in  the  model  results.  E)  Porosity  increase  in  the  brittle 
layer  at  a  MAR-type  slow-spreading  ridge-transform  system.  The  red  line  with  open  triangles 
indicates  un-cooled  600°C  isotherm  and  blue  line  with  open  triangles  indicates  hydrothermal  ly 
cooled  600°C  isotherm.  Black  dashed  line  indicates  the  effect  of  increasing  porosity  only  to  6  km 
depth.  F)  Scrpcntinization  of  the  mantle  at  an  EPR-type  system  (green  line  with  solid  green 
squares).  100%  and  25%  scrpcntinization  is  shown,  and  the  shaded  region  marks  the  serpentine 
stability  field  for  the  hydrothcrmally  cooled  ridge-transform  system.  The  model  illustration 
demonstrates  that  extent  of  serpentine  stability,  the  500°C  isotherm,  for  both  an  un-cooled  (red 
dashed  line)  and  hydrothermally  cooled  (blue  dashed  line)  thermal  model.  The  depth  of  the 
500°C  isotherm  represents  the  average  depth  calculated  by  the  3D  thermal  model  within  the 
Siqueiros  transform.  G)  Scrpcntinization  of  the  mantle  at  an  MAR-type  system  (green  line  with 
open  squares).  The  kinks  visible  in  the  100%  and  to  a  lesser  extent  the  25%  curve  result  when  the 
thickness  of  the  serpentine  layer  exceeds  6  km  depth,  the  replacement  thickness  of  layer  3. 
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Supplementary  Figure  8.  A  schematic  of  the  calculation  of  ARMBAt.r.  ARMBAt.r  = 
[(RMBAr-  RMBAr1)  +  (RMBA,  -  RMBAR2)j/2,  where  RMBAr,  and  RMBAR2  are  the  average 
RMBA  values  on  the  adjacent  ridge  segments  1  and  2,  respectively,  and  RMBAT  is  the  average 
RMBA  value  along  the  transform  fault  domain.  The  ridge  RMBA  values,  RMBAR!  and 
RMBAri,  were  calculated  by  averaging  the  values  at  the  centre  of  each  adjacent  ridge  segment 
within  a  20  km-wide  and  50  km-long  box  shown  by  green  boxes.  At  the  East  Pacific  Rise,  Juan 
de  Fuca  Ridge,  and  Southeast  Indian  Ridge,  where  the  transform  fault  zones  are  wider  due  to 
segmentation,  RMBAT  was  calculated  by  averaging  over  a  20  km-wide  swath  centred  20  km  from 
the  RTls,  indicated  by  the  blue  boxes.  At  the  Mid-Atlantic  Ridge  and  Southwest  Indian  Ridge  the 
RMBA  values  were  averaged  over  a  10  km-wide  swath  centred  20  km  from  the  RTIs.  This 
convention  is  also  used  in  determining  ARMBA,.r  for  the  2D  models. 
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S2.  Supplementary'  Tables 


Supplementary  Table  1.  Transform  fault  characteristics6'1 12 ' 24 


Transform 

fault 

Ridge 

System 

Appx. 

Latitude 

Full-slip  rate 
(cm/yr) 

Length 

(km) 

No.  of 
Segments 

ARM  BAj.r 
(mGals) 

Garrett* 

EPR 

12°S 

14.9 

130 

4 

-6.0 

Wilkes 

EPR 

8°S 

14.5 

100 

3 

-7.7 

Siqueiros 

EPR 

8°N 

11.8 

150 

5 

-13.6 

Clipperton 

EPR 

1 0°N 

11.5 

90 

1 

-3.6 

SEIR  1 

SEIR 

48°S 

7.54 

135 

1 

-8.2 

SEIR2 

SEIR 

50°S 

7.44 

80 

l 

-13.2 

Blanco 

JdFR 

42°N 

5.8 

350 

5 

-39.7 

Kane 

MAR 

24°N 

2.5 

160 

1 

6.6 

Atlantis 

MAR 

30°N 

2.4 

70 

1 

35.4 

Hayes 

MAR 

34°N 

2.2 

90 

1 

13.7 

*  Asymmetry  of  slip  at  Garrett  transform 
fault. 

^  has  been  incorporated  in  the 

thermal  model  of  Garrett  transform 

Supplemental  Table  2.  Parameters  used  in  RMBA  calculation. 


Parameter _ Value _ 

Sea  water  density,  kg  m  1 ,030 

Reference  crustal  density,  kg  nv  2,730 

Reference  mantle  density,  kg  nv  3,300 

Serpentine  density  (100%  serpentinization),  kg  nv  2,550 

Reference  crustal  thickness,  km  6 

Thermal  boundary  layer  thickness,  km  100 

Asthenospheric  temperature,  °C  1,350 

Thermal  diffusivity,  m2s'1  10'6 

Coefficient  of  thermal  expansion,  °C1  2.4  x  10'5 


59 


S3.  Supplementary  References 


1.  Blackman,  D.  K.  &  Forsyth,  D.  W.  Isostatic  compensation  of  tectonic  features  of  the 
Mid- Atlantic  Ridge  -  25°-27°30’S.  J.  Geophys.  Res.  96,  1 1741-1 1758  (1991). 

2.  Bruguier,  N.  J.,  Minshull,  T.  A.  &  Brozena,  J.  M.  Morphology  and  tectonics  of  the 
Mid-Atlantic  Ridge,  7°-12°S.  J.  Geophys.  Res.  108  (2003). 

3.  Fujiwara,  T.  et  al.  Crustal  evolution  of  the  Mid-Atlantic  Ridge  near  the  Fifteen- 
Twenty  Fracture  Zone  in  the  last  5  Ma.  Geochem.  Geophys.  Geosyst.  4  (2003). 

4.  Georgen,  J.  E.,  Lin,  J.  &  Dick,  H.  J.  B.  Evidence  from  gravity  anomalies  for 
interactions  of  the  Marion  and  Bouvet  hotspots  with  the  Southwest  Indian  Ridge: 
Effects  of  transform  offsets.  Earth  Planet.  Sci.  Lett.  187,  283-300  (2001). 

5.  Muller,  M.  R.,  Minshull,  T.  A.  &  White,  R.  S.  Crustal  structure  of  the  Southwest 
Indian  Ridge  at  the  Atlantis  II  Fracture  Zone.  J.  Geophys.  Res.  105,  25809-25828 
(2000). 

6.  Goff,  J.  A.,  Fornari,  D.  J.,  Cochran,  J.  R.,  Keeley,  C.  &  Malinverno,  A.  Wilkes 
transform  system  and  nannoplate.  Geology  21,  623-626  (1993). 

7.  Hekinian,  R.,  Bideau,  D.,  Cannat,  M.,  Francheteau,  J.  &  Hebert,  R.  Volcanic  activity 
and  crust  mantle  exposure  in  the  ultrafast  Garrett  transform-fault  near  1 3°28'S  in  the 
Pacific.  Earth  Planet.  Sci.  Lett.  108,  259-275  (1992). 

8.  Fornari,  D.  J.  et  al.  Structure  and  topography  of  the  Siqueiros  transform-fault  system  - 
Evidence  for  the  development  of  intra-transform  spreading  centers.  Mar.  Geophys. 
Res.  11,263-299(1989). 

9.  Cochran,  J.  R.,  Fornari,  D.  J.,  Coakley,  B.  J.,  Herr,  R.  &  Tivey,  M.  A.  Continuous 
near-bottom  gravity  measurements  made  with  a  BGM-3  gravimeter  in  DSV  Alvin  on 
the  East  Pacific  Rise  crest  near  9  degrees  31  '  N  and  9  degrees  50  '  N.  104,  10841- 
10861  (1999). 

10.  West,  B.  P.  &  Sempere,  J.  C.  Gravity  anomalies,  flexure  of  axial  lithosphere,  and 
along-axis  asthenospheric  flow  beneath  the  southeast  Indian  ridge.  Earth  Planet.  Sci. 
Lett.  156,  253-266(1998). 


60 


1 1.  Dziak,  R.  P.  et  al.  Recent  tectonics  of  the  Blanco  Ridge,  eastern  bianco  transform 
fault  zone.  Mar.  Geophys.  Res.  21,  423-450  (2000). 

12.  Bonatti,  E.  Vertical  tectonism  in  oceanic  fracture  zones.  Earth  Planet.  Set.  Lett.  37, 
369-379  (1978). 

13.  Karson,  J.  A.,  Tivey,  M.  A.  &  Delaney,  J.  R.  Internal  structure  of  uppermost  oceanic 
crust  along  the  Western  Blanco  Transform  Scarp:  Implications  for  subaxial  accretion 
and  deformation  at  the  Juan  de  Fuca  Ridge. ./.  Geophys.  Res.  107  (2002). 

14.  Blakely,  R.  J.  Potential  Theoty  in  Gravity  and  Magnetic  Applications  (Cambridge 
University  Press,  New  York,  NY,  1996). 

15.  Stevenson,  J.  M.,  Hildebrand,  J.  A.,  Zumberge,  M.  A.  &  Fox,  C.  G.  An  ocean-bottom 
gravity  study  of  the  southern  Juan-De-Fuca  Ridge.  J.  Geophys.  Res.  99,  4875-4888 
(1994). 

16.  Johnson,  H.  P.,  Pruis,  M.  J.,  Van  Patten,  D.  &  Tivey,  M.  A.  Density  and  porosity  of 
the  upper  oceanic  crust  from  seafloor  gravity  measurements.  Geophys.  Res.  Lett.  27, 
1053-1056  (2000). 

17.  Escartin,  J.,  Hirth,  G.  &  Evans,  B.  Nondilatant  brittle  deformation  of  serpentinites: 
Implications  for  Mohr-Coulomb  theory  and  the  strength  of  faults.  J.  Geophys.  Res. 
102,  2897-2914(1997). 

18.  Behn,  M.  D.  &  Kelemen,  P.  B.  Relationship  between  seismic  P-wave  velocity  and  the 
composition  of  anhydrous  igneous  and  meta-igneous  rocks.  Geochem.  Geophys. 
Geosyst.  4  (2003). 

19.  Ulmer,  P.  &  Trommsdorff,  V.  Serpentine  stability  to  mantle  depths  and  subduction- 
related  magmatism.  Science  268,  858-861  (1995). 

20.  deMartin,  B  ,  Hirth,  G.  &  Evans,  B.  Experimental  constraints  on  thermal  cracking  of 
peridotite  at  oceanic  spreading  centers  in  Mid-Ocean  Ridges:  Hydrothermal 
Interactions  Between  Lithosphere  and  Oceans  (eds.  German,  C.  R.,  Lin,  J.  &  Parsons, 
L.  M.)  318  (American  Geophysical  Union,  Washington,  DC,  2004). 

2 1 .  Abercrombie,  R.  E.  &  Ekstrom,  G.  Earthquake  slip  on  oceanic  transform  faults. 
Nature  410,  74-77  (2001). 


61 


22.  Canales,  J.  P.,  Detrick,  R.  S.,  Lin,  J.,  Collins,  J.  A.  &  Toomey,  D.  R.  Crustal  and 
upper  mantle  seismic  structure  beneath  the  rift  mountains  and  across  a  nontransform 
offset  at  the  Mid-Atlantic  Ridge  (35  degrees  N).  J.  Geophys.  Res.  105,  2699-2719 
(2000). 

23.  Conder,  J.  A.,  Forsyth,  D.  &  Parmentier,  E.  M.  Asthenospheric  flow  and  asymetry  of 

the  East  Pacific  Rise,  MELT  area.  ./.  Geophys.  Res.  107,  2344, 

doi:  10.1 029/200 1JB000807  (2002). 

24.  Lin,  J.,  Purdy,  G.  M.,  Schouten,  H.,  Sempere,  J.-C.  &  Zervas,  C.  Evidence  from 
gravity  data  for  focused  magmatic  accretion  along  the  Mid-Atlantic  Ridge.  Nature 
344,  627-632  (1990). 


62 


Chapter  3:  The  effects  of  mantle  rheology  and  fault  segmentation  on  melt 
generation  and  extraction  beneath  oceanic  transform  faults 


Abstract 

We  examine  mantle  melting  and  melt  extraction  beneath  fast-slipping,  segmented 
oceanic  transform  fault  systems.  Three-dimensional  (3-D)  mantle  flow  and  thermal 
structure  are  calculated  using  a  temperature  dependent  rheology  that  incorporates  a  visco¬ 
plastic  approximation  for  the  brittle  lithosphere.  Thermal  solutions  are  combined  with 
the  near-fractional,  polybaric  melting  model  of  Kinzler  and  Grove  [  1992a,  b;  1993]  to 
determine  extents  of  melting  as  well  as  major  element  melt  geochemistry.  We  investigate 
the  mantle  source  region  of  intra-transform  spreading  centers  using  the  melt  migration 
model  of  Sparks  and  Parmentier  [1991]  for  two  end  member  pooling  models:  (1)  a  wide 
mantle  region  that  incorporates  all  of  the  melt  focused  to  the  intra-transform  spreading 
center;  and  (2)  a  narrow  pooling  region  that  assumes  melt  will  not  migrate  across  a 
transform  fault  or  fracture  zone.  Models  of  a  100  km-long  transform  fault  bisected  by  a 
single  intra-transform  spreading  center  (TTSC)  reveal  that  crustal  thickness  at  the  ITSC  is 
primarily  dependent  on  the  pooling  model  and  the  mantle  potential  temperature,  Tp. 
Assuming  that  melt  is  pooled  into  the  transform  fault  domain  over  a  large  region  of  the 
surrounding  mantle  (model  1),  we  find  that  crustal  production  is  enhanced  at  ITSCs  by  up 
to  3-6.5  km.  These  model  predictions  may  explain  the  systematic  crustal  thickness 
excesses  observed  at  intermediate  and  fast-slipping  transform  faults  [Gregg  et  at .,  2007], 
Furthermore,  we  find  that  the  variability  in  major  element  lava  composition  is  primarily 
controlled  by  Tp.  Applying  these  techniques  to  the  Siqueiros  transform  we  find  that  both 
the  visco-plastic  rheology  and  a  wide  melt  pooling  region  are  required  to  explain  the 
observed  variations  in  both  crustal  thickness  and  major  element  geochemistry.  Finally, 
we  find  that  Tp  =  1350°C  and  fractional  crystallization  at  depths  of  9-15.5  km  fit  the 
majority  of  the  geochemical  observations  at  the  EPR  9°N  segment  and  the  Siqueiros 
transform  fault  system. 


To  be  submitted  as  P.  M.  Ciregg,  M.  D.  Behn,  J.  Lin,  T.  L.  Grove,  The  effects  of  mantle  rheology  and  fault 
segmentation  on  melt  generation  and  extraction  beneath  oceanic  transform  faults,  J.  Geophys.  Res. 
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1.  Introduction 


Transform  faults  are  a  first  order  feature  of  the  global  mid-ocean  ridge  system. 
Transform  fault  offsets  range  from  50  km  to  >  350  km  in  length  and  greatly  influence  the 
thermal  structure,  mantle  flow  dynamics,  and  melting  at  mid-ocean  ridges  [Fox  and  Gallo , 
1984;  Par  mender  and  Forsyth ,  1985;  Phipps  Morgan  and  Forsyth,  1988],  Investigations 
of  transform  faults  suggest  that  cool,  deep  lithosphere  at  transform  fault  offsets 
effectively  divides  mid-ocean  ridges  into  unique  segments  focusing  melt  toward  segment 
centers  and  prohibiting  shallow  across- transform  mantle  flow.  Gravity  and  seismic 
studies  at  slow-spreading  ridges  indicate  along-segment  variations  in  crustal  thickness 
with  enhanced  magma  accretion  at  ridge  segment  centers  and  crustal  thinning  towards 
transform  and  non-transform  offsets  [ Kuo  and  Forsyth,  1988;  Lin  et  al.,  1990;  Lin  and 
Phipps  Morgan,  1992;  Tolstoy  et  al.,  1993;  Detrick  et  al.,  1995;  Escartin  and  Lin,  1995; 
Canales  et  al.,  2000;  Hoqft  et  al.,  2000],  Juxtaposing  the  cold,  thick  lithosphere  of  a 
transform  fault  next  to  a  warm  steady-state  ridge  segment  also  influences  the  resultant 
lava  chemistry.  The  geochemical  transform  fault  edge  effect  has  been  noted  by  several 
studies  and  suggests  that  the  cooler  mantle  region  surrounding  the  transform  fault  system 
promotes  deeper  and  lower  extents  of  melting  near  segment  ends  [Bender  et  al,  1978; 
Bender  et  al.,  1984;  Langmuir  and  Bender,  1984;  Langmuir  et  al.,  1986;  Reynolds  et  al, 
1992;  Reynolds  and  Langmuir,  1997], 

Recent  analyses  of  residual  mantle  Bouguer  gravity  anomalies  (RMBA)  along 
several  intermediate  and  fast-slipping  transform  fault  systems  indicate  that  for  spreading 
rates  >  ~  50  mm/yr  full-rate,  the  observed  RMBA  values  within  the  transform  fault 
domain  are  more  negative  than  their  adjacent  ridge  segments  [Gregg  et  al,  2007],  The 
more  negative  RMBA  calculated  within  the  transform  fault  domain  at  fast-slipping 
transform  fault  systems  indicates  a  mass  deficit  within  these  settings.  If  the  mass  deficit 
is  due  to  crustal  thickness  variations,  these  findings  contradict  the  classic  model  for 
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transform  faults  as  regions  of  thin  crust  relative  to  adjacent  ridge  segments.  One 
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Figure  1.  Comparison  of  bathymetry  and  crustal  thickness  predictions  for  Siqueiros 
transform  fault  along  the  EPR.  (a)  EPR  segmentation  and  location  of  transform  faults. 
Note  that  all  of  the  transform  faults,  except  for  Clipperton,  are  segmented  by  at  least  one 
1TSC.  (b)  Siqueiros  transform  fault  bathymetry  with  structural  interpretation  [Fornari  el 
at.,  1989],  Solid  black  lines  indicate  the  seafloor  fabric,  circles  show  locations  of 
seamounts,  and  dashed  black  lines  show  the  locations  of  the  fracture  zones  and  transform 
fault,  (c)  Lateral  variations  in  crustal  thickness  derived  from  RMBA  calculations. 
Excess  crust  is  defined  as  the  deviation  for  the  reference  crustal  thickness  of  6  km, 
assuming  seawater,  crust  and  mantle  densities  of  1030,  2730,  and  3300  kg  nr, 
respectively. 
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explanation  for  these  observations  is  that  long-lived  spreading  centers  within  a  transform 
fault  system  may  focus  mantle  upwelling  beneath  the  transform  domain,  warm  the  fault 
system  and  result  in  excess  crustal  accretion. 

Oceanic  transform  faults  along  fast-spreading  ridges  are  often  characterized  by 
segmentation  by  intra-transform  spreading  centers  [Menard,  1967;  Menard  and  Atwater, 
1969;  Searle,  1983;  Fox  and  Gallo,  1984],  In  the  case  of  several  transform  fault  systems 
along  the  fast-spreading  East  Pacific  Rise  (EPR)  (e  g.,  Siqueiros,  Garrett,  Gofar, 
Discovery,  and  Quebrada  transform  systems),  fresh  basalts  sampled  on  ITSCs  and 
detailed  magnetic  data  suggest  accretion  is  actively  occurring  [Fornari  et  ah,  1989; 
Car  hot  te  and  Macdonald,  1992;  Hekinian  et  ai,  1992;  Hekinian  et  al,  1995;  Perfit  et  al, 
1996;  Nagle  et  al.,  2007],  Segmentation  of  oceanic  transform  faults  is  generally  thought 
to  be  the  result  of  transtensional  forces  imposed  on  the  transform  fault  zone  by  plate 
motion  reorganization  of  the  mid-ocean  ridge  system  [Menard  and  Fisher,  1958;  Menard 
and  Atwater,  1969],  Specifically,  the  change  in  spreading  direction  enables  mantle 
upwelling  beneath  the  transform  fault  domain  and  the  formation  of  small  intra-transform 
spreading  centers  (ITSC)  [Searle,  1983;  Fornari  et  al.,  1989;  Pockalny,  1997;  Pockalny 
et  al.,  1997], 

In  this  study,  we  develop  3-D  models  of  ridge-transform-ridge  systems  to  explore 
the  effect  of  transform  segmentation  on  mantle  flow,  melt  generation,  and  melt 
extraction.  Recent  modeling  investigations  of  ridge-transform-ridge  systems  have 
illustrated  the  importance  of  incorporating  visco-plastic  rheology  into  thermal  models  of 
these  settings  [Behn  et  al.,  2007],  We  compare  the  results  of  models  with  both  a  constant 
viscosity  rheology  and  a  visco-plastic  rheology  that  approximates  a  brittle  lithosphere 
[Chen  and  Morgan,  1990;  Behn  et  ai,  2007],  Calculated  thermal  structures  are  used  in 
conjunction  with  the  fractional  melting  model  of  Kinzler  and  Grove  [1992a,  b;  1993]  and 
fractional  crystallization  model  of  Yang  et  al.  [1996]  to  calculate  crustal  thickness  and 
predict  major  element  lava  compositions.  First,  crustal  thickness  and  lava  composition 
are  calculated  in  2-D  and  compared  to  observations  at  9°N  on  the  EPR  to  calibrate 
parameters  such  as  mantle  potential  temperature  and  the  maximum  distance  of  melt 
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migration.  Second,  we  investigate  a  suite  of  generic  3-D  transform  fault  models 
segmented  by  a  single  ITSC.  Finally,  we  consider  a  case  appropriate  for  the  Siqueiros 
transform  fault  system  (Figures  lb,  c)  and  compare  our  model-predicted  crustal 
thicknesses  and  lava  composition  to  geophysical  and  geochemical  observations  in  this 
region. 

2.  Model  setup 

The  3-D  mantle  flow  field  and  temperature  are  calculated  using  the  finite  element 
modeling  software  COMSOL  3.3  Multiphysics.  Mantle  flow  is  driven  by  the  divergence 
of  two  plates  each  moving  with  a  half-spreading  rate,  U0  (Figure  2).  A  stress  free 
boundary  is  assumed  at  the  base  of  the  model  (100  km  depth)  allowing  for  convective  flux 
from  the  mantle  below.  The  sides  of  the  models  are  also  stress  free.  Potential 
temperature  at  the  surface  (7'c.)  and  the  base  ( Tp )  of  the  model  are  set  to  0°C  and  Tp 
respectively.  The  model  geometry  is  varied  to  explore  the  effect  of  fault  segmentation  by 
imposing  a  single  ITSC  of  length,  L/tsc,  which  bisects  a  transform  fault  of  length  l.r 
(Figure  2).  All  model  parameters  are  given  in  Table  1 . 

2.  /  Numerical  Approach 

Conservation  of  mass,  momentum,  and  energy  are  given  by: 


V  •  u  =  0 

(1) 

P(?  +  U ' Vu)  =  "Vp  +  ri'"'S/2u 

(2) 

pC  pu  •NT  +  kS/'T  =  0 

(3) 

where  p  is  the  density  of  the  mantle,  p  is  pressure,  u  is  the  velocity  field,  i]lo,  is  the  mantle 
viscosity,  Cp  is  the  heat  capacity,  T  is  temperature,  and  k  is  thermal  conductivity. 
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Figure  2.  3D  setup  for  modeling  mantle  temperature  structure  beneath  a  segmented 
transform  fault.  Spreading  velocities  are  imposed  at  the  surface.  Other  model  parameters 
include  temperatures  at  the  surface,  Tc,  and  at  the  base  of  the  model  box,  Tp\  the  length  of 
transform,  Lj,  and  the  length  of  intra-transform  spreading  center,  Lrrsc- 
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Mantle  flow  and  temperature  are  both  largely  dependent  on  the  assumed  mantle 
viscosity.  To  explore  the  importance  of  incorporating  a  visco-plastic  rheology,  we 
compare  isoviscous  solutions  to  those  calculated  using  a  temperature-dependent  viscosity 
that  incorporates  a  visco-plastic  approximation  for  brittle  failure  [Chen  and  Morgan , 
1990;  Behn  et  a/.,  2007],  The  effective  mantle  viscosity  is  defined  by: 


Vhh 


\  But  1 lltyr  >70 ) 


(4) 


where  \},d  is  the  temperature  dependent  viscosity,  is  a  brittle  strength  approximation 
using  Byerlee’s  law,  and  )]o  is  the  maximum  viscosity  (1023  Pa-s).  The  temperature 
dependence  of  viscosity  is  approximated  by: 


/ 


nld  =  ?70exP 


\\ 


R\T  T. 


(5) 


where  Q  is  the  activation  energy  and  R  is  the  universal  gas  constant.  The  viscosity 
associated  with  brittle  failure  is  approximated  by  [Chen  and  Morgan ,  1990J: 


»7  hyr 


(6) 


where  e„  is  the  second  invariant  of  the  strain  rate  tensor,  and  the  maximum  shear  stress, 
S,  is  defined  by: 

S  =  |  {Pr~Pj)> 

where  p,  is  the  lithostatic  pressure,  and  p»  is  the  hydrostatic  pressure. 

The  incorporation  of  a  visco-plastic  rheology  promotes  focused  mantle  upwelling 
in  areas  with  high  strain  rates  [Behn  el  al.,  2007],  This  effect  is  most  pronounced  beneath 
the  ridge  axis  and  beneath  the  transform  fault  domain  where  strain  rates  are  particularly 
high.  Therefore,  models  containing  a  visco-plastic  rheology  result  in  a  warmer  thermal 
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structure  beneath  the  ridge  axis  and  transform  fault  domain  than  those  calculated  with  a 
constant  viscosity  rheology  (Figure  3a). 

Hydrothermal  circulation  within  the  upper  crust  will  promote  more  effective  heat 
loss  and  therefore  decrease  temperatures  in  the  upper  mantle  [ Phipps  Morgan  el  al. , 
1987;  Phipps  Morgan  and  Chen ,  1993;  Baker  et  al,  1996;  Maclennan  et  al,  2005].  This 
effect  is  likely  to  be  particularly  important  at  oceanic  transform  faults  where  extensive 
brittle  deformation  increases  permeability  and  produces  pathways  for  fluids  to  reach 
greater  depths.  We  assume  that  hydrothermal  circulation  is  present  at  temperatures  < 
600°C  and  depths  <  6  km  [Phipps  Morgan  et  al,  1987;  Phipps  Morgan  and  Chen ,  1993], 
Where  hydrothermal  circulation  is  assumed  to  occur  we  increase  the  thermal  conductivity 
by  an  enhancement  factor  or  Nusselt  number,  Nu  [Phipps  Morgan  et  al. ,  1987],  Because 
permeability  most  likely  decreases  with  depth  due  to  increased  confining  pressure  and 
temperature,  the  enhancement  in  conductivity  is  assumed  to  decay  exponentially  as  the 
depth  and  temperature  approach  the  cutoff  values  for  fluid  circulation: 


k  =  k0  +  k0  *  (Nu  - 1 )  *  exp 

/  T  V 

A*  h> 

*exp 

/ 

A* 

('-/ ) 

\  \  *  max  /, 

\  *  max  /  / 

where  k  is  the  effective  thermal  conductivity,  k0  is  the  reference  thermal  conductivity,  A 
is  a  smoothing  factor,  Tmax  is  the  cutoff  maximum  temperature  (600°C),  z  is  depth,  and 
Zmax  is  the  cutoff  maximum  depth  (6  km). 

2.2  Melt  Generation  and  Extraction 

Melting  within  the  mantle  at  mid-ocean  ridges  is  dependent  on  pressure  and 
temperature  as  well  as  the  chemical  composition  of  the  mantle  [Grove  et  al,  1992; 
Langmuir  et  al,  1992;  Asimow  and  Langmuir ,  2003],  As  such,  we  utilize  the  fractional 
melting  model  of  Kinzler  and  Grove  [1992b;  a;  1993],  in  which  mantle  melting  is  a 
function  of  temperature,  pressure,  and  mantle  composition.  The  benefit  of  this  method  is 
that  the  model  simultaneously  tracks  residual  mantle  and  melt  compositions,  as  well  as 
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Figure  3.  Comparison  of  melting  models  for  two  mantle  rheologies,  (a)  Comparison  of 
mantle  temperature  sampled  beneath  the  oceanic  spreading  center  for  constant  and  visco¬ 
plastic  rheology  models.  Calculated  melt  fraction  for  constant  viscosity  rheology  model 
(b)  and  visco-plastic  rheology  model  (c).  The  effect  of  latent  heat  of  melting  has  been 
incorporated  into  all  melting  models.  KG92  -  Kinzler  and  Grove  [1992a,  b;  1993],  the 
solid  line;  RJ81  Reid  and  Jackson  [1981],  the  dashed  line;  and  K03  Katz  et  at.  [2003], 
the  dotted  line.  Uo  =  50  mm/yr,  Tp  =  1 350°C,  Nu  =  4,  and  latent  heat  of  melting,  L  =  400 
J/kg. 
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extents  of  melting.  The  extent  of  melting  and  resultant  lava  composition  are  predicted 
using  a  3-step  process:  (1)  melt  production  and  composition  are  calculated  at  each  point 
in  the  melting  region  using  polybaric  fractional  melting;  (2)  melt  migrates  to  a  point  at 
the  top  of  the  melt  surface  and  pools;  and  (3)  the  pooled  melt  undergoes  fractional 
crystallization  in  the  shallow  mantle  lithosphere  and  lava  composition  is  predicted. 

2.2.1  Step  1:  Fractional  Melting 

The  polybaric  fractional  melting  model  of  Kinzler  and  Grove  (hereafter  “KG92”) 
utilizes  a  parameterized  solidus  derived  from  experimental  data  to  calculate  extents  of 
melting  [Kinzler  and  Grove,  1992b;  a;  1993],  Specifically,  in  the  spinel  stability  field 
(pressure,  p  >  9  kbars),  the  solidus  temperature,  Ts,  is  given  by  the  following 
parameterization: 

Ts=  1157+  \6*{p  -  0.001)  -  38.7*(1  -Mgtf)  -  m*(NaKtt)  -  0*(TiO2),  (9) 

and  in  the  plagioclase  stability  field  (p  <  9  kbars): 

Ts  =  1 242  +  9*(p  -  0.00 1 )  -  1 20*(  1  -  Mgtf)  -  S9*{NaKtf)  -  6.6*(Ti02).  ( 1 0) 

The  oxides  Ti02,  Na20,  CaO,  and  K20  are  estimated  with  the  nonmodal  batch 
melting  equation: 

C"'(Di<  +  F(l-P„))’  (ll) 

where  Qm  is  the  concentration  of  oxide  (i)  in  the  melt  (m),  Co  is  the  initial  mantle 
concentration,  D,:S  is  the  bulk  partition  coefficient  of  oxide  (/)  in  the  initial  solid  (5),  Pim 
is  the  bulk  partition  coefficient  of  the  oxide  in  the  melt,  and  F  is  the  melt  fraction.  The 
oxide  concentrations  are  calculated  iteratively  using  a  bisection  method.  [MgO\  and 
[FeO],  are  also  calculated  iteratively  by  solving  the  mass  balance  equations: 


(12) 


[Ft0]t-[Ff0],F-[FC0]r(l-F), 


(13) 


where  [MgO],,,  is  the  wt.  %  of  MgO  in  the  melt,  [MgO]b  is  the  wt.  %  MgO  of  the  initial 
mantle  composition,  and  [MgO],  is  the  wt.  %  MgO  of  the  residual  mantle.  Similarly, 
[FeO]m  is  the  wt.  %  of  FeO  in  the  melt,  [FeO]b  is  the  wt.  %  FeO  of  the  initial  mantle 
composition,  and  [FeO],-  is  the  wt.  %  FeO  of  the  residual  mantle.  Once  the  oxides 
concentrations  are  determined  using  equations  ( 1 1 )  -  ( 1 3),  MgHn,  and  NaKlim,  the  Mg  and 
NaK  numbers  of  the  melt,  are  calculated  using  the  following  equations: 


Mg*m- 


[MgO]m 

[MgO\n  +  [FeO]m 


(14) 


NaK#m 


(Na20+K20) 

(. Na20  +  K20  +  Ca0 ) 


(15) 


It  is  widely  accepted  that  the  process  of  melt  generation  beneath  a  mid-ocean 
ridge  spreading  center  is  adiabatic  and  occurs  over  a  range  of  pressures  [Klein  and 
Langmuir ,  1987],  Therefore,  we  utilize  the  polybaric  incremental  batch  melting  approach 
with  incomplete  melt  withdrawal  (Model  3  in  KG92)  to  approximate  fractional  melting  of 
the  mantle.  This  approach  assumes  a  constant  melt  production  rate  of  1%  melt  per  1  kbar 
of  ascent.  At  each  melting  increment,  99%  of  the  generated  melt  is  removed  and 
accumulated  elsewhere;  the  remaining  melt  and  depleted  mantle  ascend  another  1  kbar 
where  they  are  further  melted  by  1%.  As  melting  proceeds,  the  latent  heat  effect  of  the 
melting  reaction  acts  to  decrease  the  mantle  temperature.  We  incorporate  the  change  in 
temperature,  A T,  at  each  step  in  the  calculation  as  a  function  of  melt  fraction,  the  latent 
heat  of  melting,  L,  and  heat  capacity,  Cp,  using  the  following  formulation  [Turcot  te, 
1982]: 

FL 

A  T-  —  .  (16) 

C  P 


73 


To  determine  the  temperature  structure  input  into  the  KG92  melting  model,  the 
mantle  potential  temperature,  Tp,  at  the  base  of  the  model  domain  is  set  to  a  value 
between  1300°-1400°C.  Note  that,  Tp  is  defined  as  the  mantle  temperature  at  the  surface 
extracted  along  an  adiabatic  gradient  of  0.5°C/km  to  the  base  of  the  model  space  at  100 
km  depth.  For  example,  a  mantle  potential  temperature  of  1350°C  corresponds  to  a 
mantle  temperature  of  1400°C  at  100  km  depth. 

In  the  original  KG92  fractional  melting  model,  pressure  controls  the  cessation  of 
melting,  and  fractional  melting  is  assumed  to  stop  at  4  kbars.  In  our  implementation  of 
KG92,  melting  begins  when  the  adiabat  crosses  the  solidus  and  continues  until 
conductive  cooling  from  the  surface  depresses  the  temperature  below  the  solidus.  We 
assume  that  the  mantle  is  described  by  a  four-phase  assemblage  (olivine,  clinopyroxene, 
orthopyroxene,  and  an  aluminus  phase,  either  spinel  or  plagioclase  depending  on  the 
pressure),  and  utilize  the  depleted  upper  mantle  composition  (DMM)  of  Workman  and 
Hart  [2005]  as  the  initial  mantle  composition. 

2.2.2  Step  2:  Melt  Migration 

To  calculate  the  final  melt  composition  and  the  resultant  crustal  production,  our 
model  must  determine  how  melt  migrates  from  the  melting  region  to  a  nearby  spreading 
center.  Previous  3-D  melt  migration  models  have  proposed  that  melt  percolating  upward 
through  the  hot,  permeable  mantle  will  pool  along  a  cooler  less  permeable  boundary  (e.g., 
the  top  of  the  melting  region,  the  base  of  the  lithosphere,  or  the  cpx-out  reaction  zone) 
and  then  migrate  laterally,  “uphill”,  along  this  boundary  to  the  ridge  axis  [Sparks  and 
Parmentier ,  1991;  Sparks  et  al. ,  1993;  Magde  and  Sparks,  1997;  Magde  et  al.,  1997],  To 
simulate  this  two-step  melt  migration  process,  we  assume  that  melt  first  migrates 
vertically  to  the  top  of  the  melting  region  and  then  migrates  laterally  up  the  slope  of  the 
melting  region  until  it  reaches  a  local  minimum  in  the  surface,  which  typically 
corresponds  to  a  nearby  spreading  center.  Melt  migration  patterns  for  the  constant  and 
visco-plastic  mantle  rheology  models  are  shown  for  a  100  km-long  transform  fault 
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Figure  4.  Depth  to  the  top  of  melting  with  chosen  pooling  areas  for  (a)  constant 
viscosity  rheology  and  (b)  visoc-plastic  rheology.  Region  1  (delineated  by  red  outline) 
pools  the  entire  area  of  melts  that  have  the  potential  to  migrate  towards  the  1TSC.  In 
contrast,  region  2  (outlined  in  blue)  assumes  that  melts  will  not  migrate  across  the 
transform  fault.  Black  arrows  indicate  the  gradient  of  the  melt  surface  and  the  direction 
of  melt  migration.  Tp  =  1 350°C,  Nu  =  4. 
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bisected  by  a  10  km-long  ITSC  (Figure  4).  The  depth  to  the  top  of  the  melt  surface  is 
contoured  and  the  gradient  of  the  surface  (black  arrows)  is  plotted  to  indicate  the  slope  of 
the  melt  surface  and  the  direction  the  melt  will  migrate. 

Uncertainties  with  the  melt  migration  model  include:  (1)  how  far  melt  will 
migrate  along  this  boundary  before  it  freezes;  and  (2)  whether  melt  will  be  extracted 
vertically  before  it  reaches  the  ITSC.  To  explore  these  processes,  we  investigate  two 
end-member  scenarios  for  melt  migration  beneath  the  transform  fault  domain.  The  first 
scenario  assumes  a  wide  pooling  region  in  which  all  of  the  melt  projected  to  migrate  to 
the  transform  domain  will  aggregate  at  the  ITSC  (solid  red  lines  in  Figure  4).  The  second 
end-member  scenario  assumes  that  melt  is  unlikely  to  migrate  across  a  transform  fault  or 
fracture  zone  due  to  the  increased  fracture  permeability  within  the  fault  zone.  In  this 
case,  melt  is  only  pooled  at  the  ITSC  from  a  narrow  region  (solid  blue  lines  in  Figure  4). 
These  models  will  be  referred  to  as  migration  Model  1  -  wide  pooling  region,  and  Model 
2  -  narrow  pooling  region  in  all  subsequent  discussions. 

Crustal  thickness  at  the  ITSC  is  calculated  by  integrating  the  melt  production  rate 
over  the  entire  pooling  area  that  is  tapped  by  the  ITSC.  The  integrated  melt  production 
rate  is  then  divided  by  the  spreading  rate  and  the  length  of  the  ITSC  resulting  in  an 
average  crustal  thickness  for  the  entire  ITSC  [ Forsyth ,  1993], 

2.2.3  Step  3:  Fractional  Crystallization 

Finally,  to  approximate  fractional  crystallization  of  the  melt  prior  to  eruption  we 
use  the  approach  of  Yang  el  al.  [1996],  This  method  assumes  isobaric  fractional 
crystallization  and  calculates  the  olivine-plagioclase-augite  fractional  crystallization  path 
as  a  function  of  melt  composition  and  pressure  between  0. 1  to  800  MPa.  In  our  models, 
crystallization  pressure  is  constrained  by  the  depth  at  which  the  mantle  no  longer  melts. 
Yang  et  al.  [1996]  uses  algebraic  relations  obtained  from  experimentally  determined 
phase  boundaries  for  olivine-clinopyroxene-quartz  and  olivine-clinopyroxene-plagioclase 
in  pseudo-ternary  projections,  and  constraints  from  mineral-melt  exchange  reactions 
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derived  from  experimentally  produced  mineral-melt  pairs.  This  model  provides  major 
element  composition  for  each  step  of  fractional  crystallization. 

The  compositional  information  calculated  with  the  Yang  et  al.  [1996]  approach  is 
used  to  create  liquid  line  of  decent  (LLD)  plots  and  to  determine  proxies  such  as  Fes.o  and 
Na8  o,  which  can  be  compared  to  observed  major  element  variations.  LLD  plots  represent 
the  crystallization  path  the  melt  will  take  as  it  cools  in  the  shallow  parts  of  the  mantle. 
Transitions  in  the  crystallizing  phase  are  illustrated  by  the  kinks  in  the  slope  of  the  lines, 
which  in  most  cases  represent  the  transitions  from  the  equilibrium  phases  of  olivine  + 
melt,  to  olivine  +  plagioclase  +  melt,  to  olivine  +  plagioglase  +  clinopyroxene  +  melt. 
LLD  predictions  for  the  oxides  FeO  and  Na;0  are  used  as  proxies  for  depth  and  extent  of 
melting  respectively. 

3.  Numerical  Modeling  Results 

3. 1  Benchmark:  Comparison  to  EPR  9°N 

Lava  compositions  are  first  calculated  in  2-D  and  compared  to  observations  of 
crustal  thickness  and  lava  composition  at  9°N  on  the  EPR  [Baliza  and  Niu ,  1992;  Canales 
et  al.,  2003]  to  calibrate  parameters  such  as  Tp  and  the  maximum  horizontal  distance  for 
melt  migration,  fVp<M/-  For  a  constant  viscosity  rheology,  crustal  thickness  calculated  for 
Tp  —  1325°C  -  1350°C  correlates  well  with  the  seismic  observations  (Figure  5a).  While 
pooling  melt  over  a  wider  region  results  in  greater  crustal  thicknesses,  incorporating 
hydrothermal  cooling  (Nu  =  4  and  8)  decreases  the  predicted  crustal  thickness  (Figure 
5a).  As  discussed  above,  the  visco-plastic  rheology  produces  a  warmer  temperature 
structure  beneath  the  ridge  axis  resulting  in  thicker  crust,  and  we  find  that  Tp  =  1 325°C 
best  fits  the  observations  at  9°N  (Figure  5b).  Incorporation  of  hydrothermal  cooling  with 
the  visco-plastic  rheology  further  focuses  mantle  upwelling  beneath  the  ridge  creating 
larger  upwelling  velocities  and  slightly  thicker  crust  (Figure  5b). 
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Figure  5.  Comparison  of  the  predicted  crustal  thickness  at  the  ridge  axis  for  various 
widths  of  the  melt  pooling  regions  (fVpoo,  =  200,  150,  100,  and  75  km)  and  hydrothermal 
cooling  (Nu  =  1,  4,  and  8)  for  (a)  constant  viscosity  rheology  and  (b)  visco-plastic 
rheology.  The  light  gray  shaded  region  indicates  the  seismic  crustal  thickness  variations 
observed  at  the  9°N  segment  of  the  EPR  [ Canales  et  al.,  2003], 
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We  also  investigate  the  sensitivity  of  Fes.o  and  Nag.o  to  Tp  and  Wpooj  (Figure  6). 
Fex.o  and  Nas  o  are  normalizations  of  wt.  %  FeO  and  wt.  %  Na20  to  8  wt.  %  MgO  [ Klein 
and  Langmuir,  1987],  Fex.o,  which  is  classically  viewed  as  a  proxy  for  the  mean  pressure 
of  melting,  increases  as  the  pressure/depth  of  melting  increases.  Since  Na  is  highly 
incompatible,  Na«  o  is  generally  used  as  a  proxy  for  the  extent  of  melting,  decreasing  as 
the  extent  of  melting  increases.  As  Tp  increases  the  depth  of  the  onset  of  melting 
consequently  increases,  resulting  in  higher  values  of  Fes.o  (Figure  6).  Additionally, 
raising  Tp  increases  the  extent  of  melting  and  thus  Naso  is  predicted  to  decrease  (Figure 
6).  However,  while  Fex.o  and  Nax.o  are  strongly  dependent  on  Tp ,  they  vary  only  slightly 
with  mantle  rheology  and  fVpoo/. 

To  calibrate  the  values  for  Tp  and  Wpoot  used  in  our  models,  we  compare  our  2-D 
ridge  results  to  the  geochemical  observations  of  Batiza  and  Niu  [  1992],  The  Batiza  and 
Nil i  [1992]  database  is  comprised  of  major  element  basaltic  glass  data  from  a 
comprehensive  dredging  program  at  -  9-10°N  along  the  EPR  and  its  flanks.  Specifically 
we  find  that,  Tp  =  1350°C,  Nu  =  4,  and  Wpoot  -  75km  results  in  a  reasonably  good  fit  to 
the  observed  FeO  and  Na20  data  (Figure  7a,  c).  The  pressure  at  the  top  of  the  melting 
triangle  varies  for  each  model  run  depending  on  Tp  (solid  black  line  in  Figure  7).  For  Tp 
=  1350°C  the  top  of  the  melting  region  is  2.4  kbars  for  the  constant  viscosity  model  and 
2.0  kbar  for  the  visco-plastic  rheology  model.  The  dotted  line  indicates  fractional 
crystallization  at  0.2  kbars. 

For  models  with  a  visco-plastic  rheology,  a  Tp  =  1325°C  provides  the  best  match 
to  the  seismically  observed  crustal  thicknesses  at  the  EPR  9°N  (Figure  5).  However,  Tp  = 
1350°C  better  predicts  the  Na20  and  FeO  values  observed  in  microprobe  analyses  of 
glass  samples  from  the  EPR  9°N.  There  is  also  a  discrepancy  in  the  prediction  of  Fes  o 
and  Nax.o.  The  Batiza  and  Niu  [1992]  data  indicate  an  Fes.o  of  -  9,  which  corresponds  to 
our  model  prediction  with  fVpuo/  =  75  km  and  Tp  =  1325°C  (Figure  6).  Furthermore,  the 
Batiza  and  Niu  [1992]  data  indicate  Naso  ~  2.4,  which  corresponds  to  our  model 
prediction  with  Tp  =  1360°C.  Some  possible  explanations  for  this  disagreement  include 
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Figure  6.  Calculated  Nag  o  and  Feg.o  as  a  function  of  mantle  potential  temperature,  Tp.  Tp 
=  1350°C  and  Nu  =  4,  for  fVpoo/  =  75  km,  100  km,  and  150  km.  (a)  Constant  viscosity 
rheology  and  (b)  visco-plastic  rheology.  Open  circles  indicate  values  from  individual 
model  runs. 
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Figure  7.  Calculated  FeO  (top  panels)  and  Na20  (bottom  panels)  vs.  MgO  for  2D 
models  of  a  ridge  with  a  half  spreading  rate  of  50  mm/yr,  Wpooi  =  75  km,  and  Nu  =  4. 
Data  points  are  from  major  element  microprobe  analyses  of  basalt  glass  samples  collected 
along  the  9°N  segment  of  the  EPR  by  Batiza  and  Niu  [1992],  (a)  Constant  viscosity 
rheology  model  with  Tp  =  1325°C.  The  solid  black  line  indicates  fractional 
crystallization  at  top  of  the  melt  column  and  the  dashed  line  indicates  crystallization  at 
0.2  kbars.  Mantle  potential  temperature,  Tp  =  1350°C.  (b)  Models  for  visco-plastic 
rheology  with  Tp  =  1325°C.  (c)  Models  for  visco-plastic  rheology  with  Tp  =  1350°C. 
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mantle  source  heterogeneities,  variations  in  the  assumed  melt  production  rate  of  1%  melt 
per  kbar,  and  differences  in  the  radius  of  melt  pooling.  Based  on  our  3-D  model  results, 
described  later,  we  believe  that  these  discrepancies  arise  from  complexities  in  the  pattern 
of  melt  migration,  which  controls  crustal  thickness  variations.  Because  Tp  =  1350°C 
provides  the  best  match  to  the  major  element  compositions  and  Tp  =  1325°C  best  matches 
the  crustal  thickness  data,  we  investigate  3-D  models  in  later  sections  with  both  of  these 
temperature  constraints. 

We  also  compare  the  melt  productivity  predicted  by  KG92  with  two  other  melt 
models,  Reid  and  Jackson  [1981]  (referred  to  as  “RJ81”)  and  Katz  et  al.  [2003] 
(hereafter  “K03”).  RJ81  incorporates  a  pressure-temperature  solidus  relationship  and 
calculates  melting  linearly  as  a  function  of  the  difference  between  the  temperature  of  the 
mantle  and  of  the  solidus  at  a  given  pressure.  K03  incorporates  the  Hirschman  [2000] 
solidus  into  a  parameterized  equation  for  mantle  melting.  KG92  differs  from  the  melting 
approximations  of  the  RJ81  and  K03  models  because  its  solidus  is  reached  at  shallower 
depths  in  the  mantle  [Reid  and  Jackson,  1981;  Katz  et  al.,  2003],  As  a  result,  the  KG92 
model  produces  lower  extents  of  melting  than  the  K03  and  RJ81  models.  Calculated 
maximum  melt  fractions,  Fmax,  for  the  KG92  model  reach  only  13.9%  for  the  constant 
viscosity  case  and  14.5%  for  visco-plastic  rheology  (Figures  3b,  c).  In  contrast,  Fmax  for 
RJ81  is  22. 1%  for  constant  viscosity  and  25.5%  for  visco-plastic  rheology,  while  Fmax  for 
K03  is  22.3%  for  constant  viscosity  and  27.4%  for  a  visco-plastic  rheology  (Figures  3b, 
c).  Mean  extents  of  melting  range  from  5-6%  for  KG92,  9-10%  for  RJ81,  and  7-9%  for 
K03.  These  differences  become  increasingly  important  when  estimating  production  of 
crust  and  the  predicted  lava  composition.  To  match  the  crustal  thicknesses  observed  at 
the  EPR  9°N  segment,  the  RJ81  and  K03  models  require  much  shorter  distances  of  melt 
migration  to  the  ridge  axis  and/or  lower  values  of  Tp  compared  to  KG92. 
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3.2  31)  Segmented  transform  fault  results 

We  investigate  the  effects  of  fault  segmentation  on  melt  production  and  migration 
by  exploring  a  series  of  models  for  different  ITSC  lengths  (Figure  2).  Using  the  ridge 
parameters  that  best  fit  the  2-D  models  at  the  EPR  9°N,  we  assume  Tp  -  1350°C,  Nu  -  4, 
and  Uo  of  50  mm/yr  for  the  constant  viscosity  models  and  Tp  =  1325°C  and  1350°C,  for 
the  visco-plastic  rheology  models.  Similar  to  the  2-D  results,  for  a  specific  value  Tp  the 
principal  factor  controlling  crustal  thickness  is  the  size  of  the  region  over  which  melt  is 
pooled  (Figure  8).  This  is  illustrated  by  investigating  the  change  in  crustal  thickness,  Atc, 
at  the  ITSC  as  compared  to  a  reference  crustal  thickness  of  6  km.  For  the  wide  melt 
pooling  (Model  1)  we  find  that  as  Litsc  is  increased,  Atc  decreases  (Figure  8).  This  is 
because,  regardless  of  their  length,  ITSCs  pool  melt  from  a  similar  sized  mantle  region. 
Hence,  when  the  pooled  melt  is  averaged  over  the  length  of  the  ITSC,  there  is  a  larger 
effect  on  predicted  crustal  thickness  for  shorter  ITSCs.  As  the  length  of  the  ITSC 
increases,  it  behaves  more  like  a  steady-state  ridge  segment.  For  narrow  pooling  (Model 
2)  as  the  ITSC  length  is  increased,  the  predicted  crustal  thickness  also  increases  (Figure 
8).  Compared  to  calculations  with  a  constant  viscosity,  the  visco-plastic  rheology  model 
promotes  focused  upwelling  in  regions  of  increased  strain  rate  beneath  the  transform  and 
thus  systematically  predicts  greater  variations  in  crustal  accretion. 

Fault  segmentation  also  influences  melt  composition  at  the  ITSCs.  For  both  wide 
and  narrow  melt  pooling  models,  calculations  with  a  constant  viscosity  predict  both  Fe^.o 
and  Naso  to  decrease  as  the  length  of  the  ITSC  increases  (Figure  9a).  These  trends 
suggest  that  the  extent  of  melting  increases  and  the  pressure  of  melting  decreases  as  the 
ITSC  length  increases.  This  is  consistent  with  the  predicted  mantle  thermal  structure, 
because  increasing  the  length  of  the  ITCS  will  eventually  create  a  warmer,  steady-state 
ridge  segment  with  higher  extents  of  melting  and  melting  to  shallower  depths. 

In  contrast,  cases  with  a  visco-plastic  rheology  predict  relatively  uniform  lava 
composition  within  the  transform  fault  domain  (Figure  9b).  This  occurs  because  the 
visco-plastic  rheology  predicts  an  overall  warmer  thermal  structure  at  the  ITSC,  which  is 
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Figure  8:  Calculated  change  in  crustal  thickness,  Atc,  as  a  function  of  Lusc,  where  Atc  is 
defined  as  the  difference  between  the  predicted  3D  crustal  thickness  at  the  ITSC  and  a 
reference  model  assuming  6  km  of  crust.  Calculations  were  done  for  (a)  constant 
viscosity  rheology  and  (b)  visco-plastic  rheology.  Model  1  incorporates  melt  from  the 
wide  pooling  region  (circle  symbols),  while  Model  2  is  for  the  narrow  melt  pooling 
region  model  (square  symbols).  Each  symbol  is  the  result  of  a  single  3D  model  run.  Top 
and  bottom  error  bars  indicate  the  range  in  the  calculated  crustal  thickness  by  using  the 
melt  surface  contour  of  0.05  and  0.1,  respectively.  Calculations  with  Tp  =  1350°C  and 
Nu  =  4  are  indicated  by  filled  symbols,  whereas  models  with  Tp  =  1325°C  and  Nu  =  4  are 
indicated  by  open  symbols. 
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Figure  9:  Calculated  Feso  and  Na8o  as  a  function  Litsc  for  (a)  constant  viscosity 
rheology  and  (b)  visco-plastic  rheology.  Model  1  is  for  wide  melt  pooling  region,  while 
Model  2  is  for  narrow  melt  pooling  region.  Gray  circles  on  solid  lines  indicate 
calculations  for  fractional  crystallization  at  3.0  kbar.  Open  circles  on  dotted  lines 
indicate  calculations  for  fractional  crystallization  at  0.2  kbar.  Each  symbol  represents  the 
result  of  a  single  model  run.  Tp  =  1 350°C  and  Nu  =  4. 
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relatively  insensitive  to  the  length  of  the  ITSC.  This  indicates  that  the  compositional 
signature  of  the  transform  fault  is  sensitive  to  mantle  rheology.  Specifically,  while 
models  with  a  constant  viscosity  rheology  predict  a  transform  fault  effect  on  ITSC  lava 
chemistry,  models  with  a  visco-plastic  rheology  do  not. 

4.  Applications  to  the  Siqueiros  Transform  Fault  System 

4. I  Tectonic  Setting 

The  Siqueiros  transform  fault  (Figure  lb)  has  been  the  site  of  several 
investigations  including  rock  dredging  and  Alvin  dives  [Perfit  et  a/.,  1996],  magnetic 
[Carbotte  and  Macdonald ,  1992]  and  gravity  [Gregg  et  al,  2007]  surveys,  and  detailed 
geologic  mapping  [Fornari  et  al. ,  1989],  This  wealth  of  data  makes  it  an  ideal  location  to 
investigate  the  effects  of  fault  segmentation.  The  Siqueiros  transform  system  is  a  -140 
km-long  left-lateral  offset  located  between  8°  20’N  and  8°  30’N  on  the  EPR.  The  half 
slip-rate  of  the  transform  system  is  -  55  mm/yr  [DeMels  et  al.,  1990;  1994],  Bathymetry 
and  side-scan  sonar  collected  from  Siqueiros  reveal  that  the  transform  system  is 
segmented  into  five  fault  strands  separated  by  four  ITSCs  labeled  A-D  from  west  to  east 
(Figure  10)  [ Fornari  et  al.,  1989],  The  current,  segmented  geometry  of  Siqueiros 
transform  is  thought  to  be  the  result  of  several  discrete  periods  of  counterclockwise 
rotation  in  the  spreading  direction  (~  4-8°)  over  the  past  3  m.y.  [Pockalny  et  al.,  1997], 

Residual  mantle  Bouguer  gravity  anomalies  (RMBA)  calculated  along  the 
Siqueiros  transform  fault  system  suggest  a  mass  deficit  within  the  transform  fault  domain 
[Gregg  et  al.,  2007],  Assuming  that  the  entire  gravity  anomaly  is  due  to  crustal  thickness 
variations  the  RMBA  data  indicate  excess  crustal  thickness  (as  compared  to  a  reference 
crust  of  6  km)  of  up  to  1.5  km  associated  with  bathymetric  high  associated  with  ITSC  D 
(Figure  2c).  Excess  crust  is  also  associated  with  ITSCs  B  and  C  along  the  central  portion 
of  the  transform  fault,  and  in  the  inactive  fracture  zones  on  either  side  of  the  Siqueiros 
transform  system  [Gregg  et  al.,  2007],  Gravity  analysis  did  not  reveal  excess  crust  at 
ITSC  A. 
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Figure  10:  Shaded  bathymetry  for  the  Siqueiros  Transform  showing  the  locations  of  rock 
dredge  and  Alvin  dive  samples  (colored  symbols)  collected  along  the  transform  fault 
[Hays,  2003],  The  solid  black  line  indicates  the  plate  boundary,  while  the  letters  denote 
individual  ITSCs. 
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We  incorporate  the  geometry  of  Siqueiros  transform  system  into  a  3-D  model  to 
evaluate  melting  and  melt  migration  beneath  the  fault  zone.  ITSCs  A-D  range  in  length 
from  10  km  at  ITSC  A  to  <  5  km  at  ITSC  D.  Due  to  the  obliquity  of  the  ITSCs  within  the 
transform  fault  and  the  computational  difficulty  of  the  exact  geometry,  the  orientations  of 
each  of  the  ITSCs  have  been  squared  off  so  that  they  are  orthogonal  to  the  transform  fault 
strands  and  parallel  to  the  adjacent  ridge  segments  (Figure  1 1). 

4.2  Crustal  Thickness  Predictions  at  Siqueiros 

Model-predicted  crustal  thickness  variations  at  the  Siqueiros  fault  system  are 
strongly  dependent  on  mantle  rheology  (Figure  11a  and  b).  The  constant  viscosity 
rheology  predicts  a  cooler  mantle  beneath  the  transform  system  and,  subsequently,  lower 
extents  of  melting  compared  to  the  visco-plastic  rheology  (Figure  11a).  Along  the 
western  portion  of  Siqueiros,  the  predicted  crustal  thickness  variations  due  to  the  wide 
melt  pooling  (Model  1)  match  the  gravity-derived  crustal  thicknesses  calculated  for 
ITSCs  A  and  B  (Figure  1  lc).  However,  the  predicted  crustal  thickness  diminishes  to  the 
east  and  does  not  match  the  observed  gravity-derived  crustal  thicknesses  at  ITSCs  C  and 
D.  Similar  to  the  single  ITSC  models,  segmentation  geometry  plays  an  important  role  in 
how  the  melt  migrates  into  the  transform  fault.  At  ITSC  A,  the  effect  of  the  western 
ridge-transform  intersection  (RTI)  results  in  a  local  minimum  in  the  top  of  the  melting 
region  that  is  offset  to  the  north  of  ITSC  A.  While  we  assume  that  these  melts  would  be 
distributed  along  ITSC  A,  the  offset  suggests  that  this  would  not  be  a  steady-state 
geometry.  Because  of  the  proximity  of  ITSC  D  to  the  eastern  RTI  and  its  relatively  short 
length  (~  5km),  ITSC  D  does  not  induce  a  local  minimum  in  the  top  of  the  melting 
region.  Therefore,  melt  produced  near  ITSC  D  is  predicted  to  migrate  to  the  eastern  RTI 
and,  thus,  will  not  contribute  to  crustal  accretion  at  ITSC  D  (Figure  1  lc). 

The  key  difference  between  the  visco-plastic  and  constant  viscosity  models  is  that 
the  visco-plastic  rheology  predicts  a  warmer  shallow  mantle  beneath  the  transform  fault 
domain  and  higher  overall  extents  of  melting.  The  warmer  temperature  structure  causes 
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Figure  11:  Results  from  the  Siqueiros  Transform  melting  models,  (a)  Depth  to  the  top 
of  melting  calculated  for  a  ridge-transform-ridge  system  with  the  geometry  of  Siqueiros 
and  a  constant  viscosity  rheology.  Arrows  indicate  the  local  spatial  derivative  of  the  melt 
surface  that  is  assumed  to  be  the  melt  migration  direction.  Red  solid  line  indicates  wide 
melt  pooling  regions  (Model  1).  Blue  solid  lines  indicate  the  narrow  melt  pooling 
regions  (Model  2),  which  assumes  that  melt  does  not  migrate  across  the  transform  fault, 
(b)  The  depth  to  the  top  of  the  melt  surface  for  the  visco-plastic  rheology  case,  (c) 
Calculated  and  estimated  crustal  thickness  values  for  the  constant  viscosity  model.  Tp  - 
1350°C  and  Nu  =  4.  Gray  bars  show  the  gravity-derived  crustal  thicknesses  at  the 
Siqueiros  transform  fault  from  Gregg  et  al.  [2007],  Red  bars  indicate  the  predicted 
crustal  thickness  for  the  wide  pooling  region  model.  Blue  bars  indicate  the  predicted 
crustal  thicknesses  for  the  narrow  pooling  region  model.  Note  that  no  crustal  accretion  is 
predicted  at  spreading  center  D,  as  it  is  not  a  local  minimum  in  the  melt  surface,  (d) 
Crustal  thickness  predictions  for  the  visco-plastic  rheology  model  Tp  =  1325°C  and  Nu 
=  4.  Note  that  this  model  does  a  much  better  job  of  predicting  the  crustal  thicknesses  as 
observed  from  the  gravity  calculations.  Crustal  thickness  predictions  for  Tp  -  1350°C  are 
indicated  by  the  dotted  lines. 
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each  ITSC  to  produce  a  local  minimum  in  the  top  of  the  melting  region  and  therefore 
enhances  crustal  accretion  at  each  ITSC  (Figure  1  lb).  For  Tp  =  1325°C  the  predicted 
crustal  thickness  for  the  wide  pooling  model  (Model  1)  generates  a  good  fit  to  the 
gravity-derived  crustal  thickness  data  at  all  4  ITSCs.  Varying  Tp  from  1325°C  to  1350°C 
(dotted  bars  in  Figure  lid)  results  in  predicted  crustal  thickness  values  that  are 
approximately  the  same  as  those  derived  from  gravity  at  ITSCs  A,  B,  and  C,  but  over¬ 
predicts  the  crustal  thickness  produced  at  ITSC  D.  Finally,  the  crustal  thickness  values 
predicted  using  narrow  melt  pooling  (Model  2)  for  both  the  constant  and  visco-plastic 
rheologies  are  significantly  smaller  than  those  derived  from  the  RMBA  data,  indicating 
that  a  wide  melt  pooling  region  is  necessary  to  produce  the  crustal  thickness  variations 
calculated  from  the  gravity  anomaly. 

4.3  Geochemical  Predictions  at  Siqueiros 

We  utilize  major  element  analyses  of  the  glassy  chill  margins  from  fresh  basalts 
collected  via  dredge  and  Alvin  dive  along  the  Siqueiros  fault  system  (Figure  10)  [ Perfit  et 
al. ,  1996;  Hays ,  2004],  Using  the  isobaric  fractional  crystallization  model  of  Yang  et  al. 
[1996],  the  predicted  liquid  lines  of  descent  (LLD)  are  compared  to  the  observed  major 
element  composition  of  the  Siqueiros  transform  lavas.  Fractional  crystallization  pressure 
is  fixed  by  the  depth  where  melting  ceases  (solid  lines  in  Figures  12  and  13)  or  at  0.2 
kbars  for  the  average  seafloor  depth  (dashed  lines  in  Figures  12  and  13). 

In  general,  the  visco-plastic  rheology  with  wide  pooling  and  Tp  =  1350°C  does  the 
best  job  of  fitting  the  majority  of  the  data  (Figures  12  and  13).  Both  the  constant  and 
visco-plastic  rheology  models  compare  favorably  to  the  observed  FeO  composition  and 
there  is  very  little  difference  between  the  wide  and  narrow  pooling  models.  The  warmer 
temperatures  resulting  from  the  visco-plastic  rheology  model  produce  lower  FeO 
implying  lower  pressures  of  mantle  melting  compared  to  the  constant  viscosity  rheology 
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Figure  12.  Liquid  lines  of  decent  (LLD)  model  calculations  of  FeO  vs.  MgO  Colored 
LLD  in  each  panel  correspond  to  the  color  of  the  data  point  for  each  ITSC  or  transform 
fault  strand  as  shown  in  the  legend  on  the  lower  right  and  in  Figure  10.  Solid  LLD  lines 
indicate  fractional  crystallization  at  the  top  of  the  melt  surface,  which  ranges  in  depth 
from  18-24  km  (6-8  kbars)  for  constant  viscosity  and  9-15.5  km  (3-5  kbars)  for  visco- 
plastic  rheology,  while  dashed  LLD  lines  indicate  fractional  crystallization  at  0.2  kbar 
which  would  be  expected  for  fractional  crystallization  in  a  shallow  magma  lens.  All 
models  were  run  for  Tp  =  1350°C  and  Nu  =  4.  The  gray  shaded  regions  indicate  the 
predicted  model  region  for  Tp  —  1325°C.  Models  were  run  for  (a)  constant  viscosity 
rheology  and  (b)  variable  rheology  models  and  for  a  wide  melt  pooling  region  (top 
panels)  and  a  narrow  melt  pooling  region  (bottom  panels). 
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Figure  13.  Liquid  lines  of  decent  (LLD)  model  calculations  of  Na^O  vs.  MgO.  Colored 
LLD  in  each  panel  correspond  to  color  of  the  data  point  for  each  ITSC  or  transform  fault 
strand  as  shown  in  the  legend  on  the  lower  right  and  in  Figure  10.  Solid  LLD  lines 
indicate  fractional  crystallization  at  the  top  of  the  melt  surface,  which  ranges  in  depth 
from  18-24  km  (6-8  kbars)  for  constant  viscosity  and  9-15.5  km  (3-5  kbars)  for  visco¬ 
plastic  rheology,  while  dashed  LLD  lines  indicate  fractional  crystallization  at  0.2  kbar 
which  would  be  expected  for  fractional  crystallization  in  a  shallow  magma  lens.  All 
models  were  run  for  Tp  =  1350°C  and  Nu  =  4.  The  gray  shaded  regions  indicate  the 
model  region  for  Tp  =  1325°C.  Models  were  run  for  (a)  constant  viscosity  and  (b)  visco¬ 
plastic  rheology  models  for  a  wide  melt  pooling  region  (top  panels)  and  a  narrow  melt 
pooling  region  (bottom  panels). 
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(Figures  12a,  b).  Fractional  crystallization  at  the  depth  of  melt  termination,  18-20  km 
and  9-15.5  km  for  the  constant  and  visco-plastic  rheologies  respectively  (solid  LLD  in 
Figures  12  and  13),  fits  the  data  much  better  than  fractional  crystallization  at  shallower 
depths  (dashed  LLD  in  Figures  12  and  13).  The  discrepancy  between  the  Siqueiros  FeO 
data  and  the  visco-plastic  rheology  predicted  LLD  may  indicate  that  crystallization  is 
polybaric  and  a  mixture  of  melts  fractionally  crystallized  at  various  depths  (Figure  12b). 
Incorporating  a  lower  Tp  of  1325°C  results  in  shallower  melting  with  lower  FeO  and 
higher  Na^O  predictions  that  do  not  agree  as  well  with  the  Siqueiros  data  for  either 
mantle  rheology.  The  constant  viscosity  model  under-predicts  the  extent  of  melting,  i.e., 
Na:0  values  that  are  too  high,  for  both  wide  and  narrow  pooling  (Figure  13a).  This  is 
consistent  with  the  under-prediction  of  crustal  thickness  using  the  constant  viscosity 
model. 

The  predicted  trends  in  Fex.o  and  Nax.o  for  the  Siqueiros  transform  system  are  very 
similar  to  those  predicted  by  the  generic  ridge-transform-ridge  models  (Figures  9a,  b). 
As  the  ITSCs  decrease  in  length  from  west  to  east  along  the  Siqueiros  transform  fault,  the 
visco-plastic  rheology  predicts  little  to  no  variation  in  the  Fex.o  and  Nax.o  (Figure  14b). 
On  the  contrary,  the  constant  viscosity  rheology  predicts  that  both  Fex.o  and  Nax.o  will 
increase  indicating  that  the  extent  of  melting  decreases  and  the  depth  of  melting 
increases.  The  observed  range  in  Nax.o  (shown  by  the  error  bars  in  Figure  14)  appears  to 
increase  slightly  along  the  transform  fault  in  agreement  with  the  model  predictions,  but 
Fex.o  decreases  slightly  from  west  to  east  along  the  transform  fault  system.  Both 
rheological  models  intersect  the  range  of  Fex.o  and  Nax.o  observed  at  all  of  the  ITSCs  and 
fault  strand  AB,  but  do  not  do  a  very  good  job  of  predicting  the  values  observed  at  fault 
strand  BC. 

5.  Discussion 

5. 1  Factors  controlling  crustal  thickness  variations 

Melt  migration  beneath  mid-ocean  ridge  spreading  centers  is  poorly  constrained. 
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Figure  14:  Fes.o  and  Nag.o  predictions  for  each  of  the  ITSCs  and  their  connecting  fault 
segments  at  the  Siqueiros  Transform.  Mean  data  values  are  indicated  by  solid  dots 
[Hays,  2004],  while  error  bars  indicate  the  total  spread  observed  in  the  data.  Dark  gray 
boxes  indicate  predictions  for  a  wide  melt  pooling  region  and  fractional  crystallization 
from  0.2  kbar  to  the  top  of  the  melting  region  (~6-8  kbar  for  constant  viscosity  and  3-5 
kbar  for  visco-plastic  rheology).  Light  gray  shaded  boxes  indicate  the  predictions  for  a 
narrow  melt  pooling  region  and  fractional  crystallization  from  0.2  kbar  to  the  top  of  the 
melting  region  ( — 6 — 8  kbar  for  constant  viscosity  and  3-5  kbar  for  visco-plastic  rheology). 
Open  boxes  show  the  predictions  for  melts  migrated  to  the  transform  fault  strands  AB 
and  BC. 
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In  2-D  ridge  models,  the  predicted  crustal  thickness  variations  are  sensitive  to  Wpooi  and 
Tp.  For  example,  for  Uo  =  50  mm/yr,  Wpooi  =  75  km,  and  Tp  =  1350°C,  the  KG92  melt 
model  predicts  a  crustal  thickness  of  6-7  km,  while  the  K03  and  RJ81  melting  models 
predict  thicker  crust  on  the  order  of  8-10  km.  Therefore,  to  match  the  crustal  thickness 
variations  observed  at  the  EPR  9°N  segment,  the  RJ81  and  K03  models  require  the 
migration  of  melt  over  shorter  distances  to  the  ridge  axis  or  lower  values  for  Tp. 

A  similar  relationship  is  observed  in  3-D  at  ITSCs  within  a  transform  fault 
system.  Previous  studies  have  assumed  that  melts  extruded  at  ITSCs  came  from  directly 
below  the  transform  fault  domain  or  from  the  nearby  ridge  axis,  and  have  discussed  the 
possibility  that  melts  may  be  drawn  in  across  the  transform  fault  from  large  regions  of  the 
surrounding  mantle  [Fornari  et  al.,  1989;  Hekinian  et  al.,  1992;  Hekinian  et  al.,  1995; 
Perfit  el  al.,  1996;  Wendt  et  al.,  1999].  If  we  assume  that  the  negative  RMBA  observed 
along  the  Siqueiros  fault  system  is  due  to  crustal  thickness  variations,  melts  must 
somehow  focus  into  the  transform  fault  from  large  regions  of  the  mantle.  The  wide  melt 
pooling  model  for  the  Siqueiros  transform  predicts  that  melts  must  travel,  on  average, 
~40  km  to  reach  an  1TSC.  This  is  approximately  half  the  distance  that  was  required  for 
melts  to  migrate  (75  km)  in  order  to  match  the  EPR  9°N  data  in  our  2-D  melting  models. 
However,  since  there  are  fresh  lava  flows  observed  along  the  fault  strands  separating  the 
ITSCs,  it  is  clear  that  some  of  the  melt  does  not  make  it  to  the  spreading  segments  [Perfit 
et  al.,  1996],  Overall,  we  conclude  that  to  generate  enhanced  crustal  thickness  within 
fast-slipping  transform  systems,  melt  must  migrate  into  the  transform  domain  from  a 
wide  region  of  underlying  mantle. 

RMBA  derived  crustal  thickness  variations  [Gregg,  et  al.,  2007]  at  the  Garrett 
transform  fault  also  suggest  that  melt  must  come  from  a  wide  region  of  the  mantle. 
However,  the  melt  generated  in  the  mantle  surrounding  the  Garrett  transform  does  not 
appear  to  migrate  directly  into  the  transform  fault  and  aggregate  at  the  ITSCs.  The 
Garrett  transform  is  -  130  km-long  with  three  ITSCs  containing  four  ridges  ~  9  km-long 
that  lie  in  between  four  fault  strands  [Hekinian  et  al.,  1992].  Peridotites  collected  from 
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the  Garrett  transform  are  very  depleted,  indicating  high  degrees  of  melting  in  the  mantle 
regions  surrounding  the  fault  system  [Niu  and  Hekinian,  1997a;  Niu  and  Hekinian, 
1997b;  Constantin ,  1999].  However,  the  majority  of  the  excess  crust  observed  at  the 
Garrett  system  is  associated  with  the  flanking  topography  correlating  to  the  transform 
fault  and  fracture  zones  and  not  with  the  ITSCs  within  the  transform  domain.  The  Garrett 
ITSCs  appear  to  correspond  to  regions  of  modest  crustal  thicknesses  equal  to  or  less  than 
the  adjacent  ridge  segments  [Gregg  et  al.,  2007],  These  observations  suggest  that  the 
melts  pooled  from  a  wide  mantle  region  surrounding  the  Garrett  system  are  extruded  on 
the  flanks  of  the  transform  fault  where  more  permeable  pathways  may  have  been 
produced  by  brittle  deformation. 

5.2  The  Siqueiros  Transform  Effect  on  Mantle  Melting 

The  geochemical  transform  fault  effect  suggests  that  the  cooler  mantle  region 
surrounding  the  transform  fault  system  will  promote  deeper  and  lower  extents  of  melting 
[Bender  et  al.,  1978;  Bender  et  al.,  1984;  Langmuir  and  Bender,  1984;  Langmuir  et  al., 
1986;  Reynolds  et  al.,  1992;  Reynolds  and  Langmuir,  1997],  However,  previous  studies 
have  illustrated  considerable  overlap  in  the  observed  composition  of  lavas  collected  along 
the  EPR  9°N  segment  and  the  Siqueiros  transform  fault  system  [Perfit  et  al.,  1996;  Hays, 
2004]  and  only  modest  along-segment  variation  in  the  composition  of  lavas  from  the  9°N 
ridge  segment  approaching  the  Siqueiros  transform  fault  [Langmuir  et  al.,  1986], 
Furthermore,  the  investigation  of  U-series  disequilibria  in  basalts  collected  at  the 
Siqueiros  transform  fault  suggests  that  the  upwelling  and  melting  rates  beneath  the 
transform  fault  are  very  similar  to  those  observed  in  samples  from  9°30’  N  on  the  EPR 
[Lundstrom  et  al.,  1999], 

The  incorporation  of  the  visco-plastic  mantle  rheology  into  the  thermal  models  for 
a  transform  fault  creates  a  warmer  thermal  structure  beneath  the  transform  fault  domain. 
This  viscosity  effect  coupled  with  a  segmented  transform  fault  geometry  results  in 
melting  processes  at  ITSCs  that  mimic  the  processes  occurring  at  nearby  ridge  segments, 
and  predicts  extents  of  melting  and  lava  composition  at  the  ITSCs  that  are  comparable  to 
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those  observed  at  a  steady-state  ridge  segment.  Therefore,  not  only  does  the  visco-plastic 
rheology  model  match  the  geophysical  and  geochemical  observations  at  Siqueiros,  this 
rheology  does  a  very  good  job  of  explaining  why  the  ITSCs  of  Siqueiros  are  producing 
EPR-like  lavas. 

5.3  High-pressure  Fractional  Crystallization 

While  high  pressure  fractional  crystallization  has  been  proposed  by  Michael  and 
Cornell  [1998]  to  explain  compositional  variability  at  slow-spreading  ridges  with  low 
magma  flux,  Langmuir  et  al.  [1992]  maintain  that  there  are  many  problems  associated 
with  this  mechanism.  The  primary  argument  of  Langmuir  et  al.  [1992]  is  that  in  a  ridge 
setting  the  lithosphere  must  be  sufficiently  thick  to  shut  melting  off  at  depth  and  that 
high-pressure  crystallization  must  be  consistent  with  melting  processes.  By  combining 
fractional  melting  and  crystallization  models  with  mantle  thermal  models  we  are  able  to 
quantify  both  the  termination  of  melting  and  the  onset  of  crystallization.  At  the  9°N 
segment  of  the  EPR  our  predicted  fractional  crystallization  depth  range  of-  0-3  km  (0.2- 
1.0  kbar)  matches  major  element  composition  observations  [Batiza  and  Niu,  1992]  and 
crustal  thickness  variations  [Canales  et  al.,  2003],  as  well  as  previous  estimations  for 
shallow  crystallization  beneath  the  EPR  ridge  axis  [e.g.,  Michael  and  Cornell.  1998].  On 
the  other  hand,  we  find  that  models  of  fractional  crystallization  at  much  greater  depths, 
9-18  km  (3-6  kbars),  fit  the  observed  variations  in  composition  at  the  Siqueiros 
transform  fault  ITSCs  [ Perfit  et  al.,  1996;  Hays,  2004],  In  fact,  models  of  shallow 
crystallization  do  a  very  poor  job  of  predicting  the  observed  major  element  variations  at 
the  Siqueiros  transform  fault.  Specifically,  if  shallow-level  fractional  crystallization  were 
occurring  we  would  expect  the  observed  compositions  to  fall  upon  the  0.2  kbar  line  in 
Figures  12  and  13.  Since  our  models  take  into  account  the  full  melting  process  and 
begin  crystallization  at  depths  above  melt  termination,  they  are  self-consistent  given  the 
physical  constraints  of  melting  and  crystallization  processes. 
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5.4  Melt  Extraction  at  Siqueiros  Transform  Fault 

The  observed  major  and  trace  element  data  from  the  Siqueiros  transform  system 
vary  widely  from  enriched  mid-ocean  ridge  basalt  (EMORB)  to  depleted  mid-ocean  ridge 
basalt  (DMORB).  However,  lavas  from  a  single  tectonic  region  within  Siqueiros  are 
relatively  homogenous  (e.g.,  Figure  13)  [Perfit  et  al.,  1996;  Hays  et  al.,  2004],  These 
findings  have  led  previous  researchers  to  suggest  that  there  are  multiple  mantle  sources 
beneath  the  transform  fault  and  that  each  ITSC  and  fault  strand  is  independent  of  its 
neighbor  and  taps  a  different  melt  lens  and/or  mantle  supply  [Fornari  et  al.,  1989;  Perfit 
et  al,  1996;  Hays,  2004],  Fornari  et  al.  [1989]  proposed  that  the  variations  in  the 
observed  lava  chemistry  result  from  melt  accumulation  at  different  depths  due  to  the 
complex  thermal  structure  of  Siqueiros.  The  range  of  normal  MORB  compositions  found 
in  the  lavas  collected  from  the  ITSCs  suggest  the  existence  of  a  frequently  replenished 
melt  lens  beneath  each  ITSC  [Perfit  et  al.,  1996].  Perfit  et  al.  [1996]  further  suggested 
that  the  primitive  lavas  collected  from  fault  strand  AB  might  be  the  result  of  tapping  into 
the  unmixed,  undifferentiated  lower  levels  of  a  magma  chamber,  possibly  from  one  of  the 
adjacent  ITSCs.  All  of  these  previous  studies  suggest  that  the  three-dimensional 
magmatic  regime  beneath  Siqueiros  is  convoluted,  and  are  unclear  as  to  how  each  ITSC 
is  supplied  from  the  mantle. 

Until  now,  the  process  of  melt  feeding  into  the  transform  fault  domain  has  been 
thought  of  as  fairly  two-dimensional  in  that  the  supply  of  melt  must  either  come  from 
directly  below  the  transform  fault  or  from  the  nearby  ridge  axis.  Our  3-D  thermal  and 
melt  models  for  Siqueiros  indicate  that  melt  present  in  the  mantle  regions  flanking  the 
transform  fault  should  also  supply  the  ITSCs.  By  using  classic  melt  migration  techniques 
[e.g..  Sparks  and  Parmentier,  1990]  that  pool  melts  from  a  large  region  of  the  mantle  and 
fractionally  crystallizing  them  at  the  depth  of  melt  termination,  we  can  reproduce  the 
same  range  of  major  element  lava  composition  that  is  observed  (Figures  12b,  and  13b). 
The  results  of  our  models  indicate  that  melts  extracted  at  Siqueiros  are  aggregated  from  a 
wide  region  of  the  mantle  and  experience  fractional  crystallization  at  depths  raging  from 
~  9-18  km.  This  suggests  that  the  variations  in  the  geochemisty  of  lavas  from  Siqueiros 
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transform  fault  can  be  explained  by  the  complexities  of  migrating  melt  into  the  transform 
fault,  and  do  not  require  (nor  preclude)  multiple  sources  beneath  Siqueiros. 

6.  Conclusions 

Our  investigation  of  melt  extraction  at  segmented  oceanic  transform  faults 
illustrates  that  melt  migration  into  the  transform  fault  domain  is  a  fully  three-dimensional 
process  that  incorporates  melts  from  a  wide  region  of  the  mantle  surrounding  the 
transform  fault.  Assuming  that  melt  is  migrated  into  the  transform  fault  domain  from  a 
large  mantle  region,  we  find  that  crustal  production  is  enhanced  at  ITSCs  by  2.5  km  for 
models  incorporating  an  isoviscous  rheology  and  up  to  3-6.5  km  for  visco-plastic 
rheology  models.  These  model  predictions  may  explain  the  systematic  excess  crustal 
thickness  observed  at  intermediate  and  fast-slipping  transform  faults  [Gregg  el  al.,  2007], 
Specifically,  to  match  the  gravity-derived  crustal  thickness  variations  observed  within  the 
Siqueiros  transform  fault  (>  1.5  km),  melt  must  aggregate  from  a  wide  pooling  area. 
Finally,  utilizing  a  visco-plastic  rheology,  with  Tp  =  I350°C,  in  conjunction  with  the 
fractional  melting  model  of  Kinzler  and  Grove  [1992a,  b;  1993]  and  the  fractional 
crystallization  model  of  Yang  et  al.  [1996],  we  demonstrate  that  the  observed  variations 
in  major  element  composition  from  the  Siqueiros  transform  fault  can  be  explained  by 
pooling  melts  from  a  wide  region  of  the  mantle  and  fractionally  crystallizing  them  at  the 
depth  of  melt  termination,  ~  9—  1 8  km. 
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9.  Tables 


Table  1.  Model  parameters 


Parameter 

Value 

A 

Smoothing  coefficient  for  hydrothermal  cooling  equation 

0.75 

cP 

Specific  heat  capacity,  J  kg-1  °K  1 

1.25  x  10 

Vo 

Reference  maximum  viscosity.  Pa  s 

1023 

k 

Thermal  diffusi vity,  m2  s'1 

10‘6 

k{) 

Reference  thermal  conductivity,  W  m  1  °K  1 

3 

L 

Latent  heat  of  melting,  J  kg  1 

400 

lt 

Transform  fault  length,  km 

100 

L/tsc 

Intra-transform  spreading  center  length,  km 

5-25 

Nu 

Nusselt  number 

1  -  8 

0 

Activation  energy,  J  mol1 

250 x  10’ 

P 

Reference  mantle  density,  kg  m  3 

3300 

R 

Universal  gas  constant,  J  mof1  °K  1 

8.3114 

Tc 

Temperature  at  the  surface  of  the  thermal  model,  °C 

0 

Tmax 

Maximum  temperature  of  hydrothermal  fluid  circulation,  °C 

600 

TP 

Mantle  potential  temperature  used  in  melt  models,  °C 

1300-  1400 

Ho 

Imposed  spreading  half-rate,  mm  yr'1 

50 

Wpoo/ 

Horizontal  length-scale  for  melt  migration,  km 

75  -  200 

Zmax 

Maximum  depth  of  hydrothermal  fluid  circulation,  km 

6 

Table  2. 

Variables 

Variable 

Co 

Initial  oxide  concentration 

c,m 

Concentration  of  oxide,  i,  in  the  melt 

D,.s 

Bulk  partition  coefficient  of  oxide,  i,  in  the  initial  solid 

Atc 

Change  in  crustal  thickness,  km 

^// 

Second  invariant  of  the  strain  rate  tensor,  s’1 

F 

Melt  fraction,  % 

F  max 

Maximum  melt  fraction,  % 

Vbyr 

Brittle  strength  viscosity,  Pa  s 

ri,d 

Temperature  dependent  viscosity.  Pa  s 

no,i 

Total  viscosity.  Pa  s 

k 

Effective  thermal  conductivity,  W  m  1  °K  1 

p 

Pressure  of  melting,  kbars 

Pr 

Lithostatic  pressure.  Pa 

Pn 

Hydrostatic  pressure,  Pa 

Pi.m 

Bulk  partition  coefficient  of  oxide,  i,  in  the  melt 

S 

Maximum  shear  stress.  Pa 

T 

Temperature,  °C 

Ts 

Solidus  temperature,  °C 

u 

Velocity  field 
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Appendix:  Seismicity  characteristics  of  the  equatorial  East  Pacific  Rise  and 
its  offsets 


Abstract 

We  utilize  hydroacoustic  data  collected  from  May  1996  to  November  2001  by 
NOAA/PMEL’s  Equatorial  Pacific  autonomous  underwater  hydrophone  (AUH)  array  to 
investigate  the  first  order  spatial  and  temporal  variability  of  seismicity  of  the  equatorial 
Pacific  within  the  array  and  in  particular,  the  East  Pacific  Rise  (EPR).  Our  statistical 
analyses  reveal  strong  evidence  of  seismic  clustering  along  the  EPR.  The  majority  of 
earthquakes  are  clustered  in  seismic  swarms  in  which  individual  events  occurred  within  a 
few  km  and  minutes  or  hours  of  each  other.  We  examine  the  seismicity  and  segmentation 
characteristics  of  eight  transform  fault  systems  in  our  study  area:  Clipperton,  Siqueiros, 
Quebrada,  Discovery,  Gofar,  Yaquina,  Wilkes,  and  Garrett  transform  faults.  While  the 
seismic  moment  release  recorded  in  the  CMT  catalog  over  the  past  30  years  increases 
with  increasing  brittle  area  for  each  transform  fault,  the  number  of  recorded  AUH  events 
does  not  appear  to  be  correlated  to  brittle  area.  On  the  other  hand,  with  the  exception  of 
Clipperton  transform  fault  and  the  G3  segment  of  Gofar  transform,  the  number  of 
recorded  AUH  events  increase  with  increasing  transform  segment  length.  More  than 
90%  of  the  seismicity  recorded  on  the  ridge  axis  is  associated  with  overlapping  spreading 
centers  (OSC).  While  many  of  these  events  can  be  attributed  to  tectonic  activity  on  the 
OSC,  there  is  also  evidence  for  magmatically  induced  seismicity. 
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Tables 


Table  1.  First-order  ridge  segmentation  within  the  study  area 


No.  of  seismic 

North  end _  South  end  events 


Segment 

Lat  °N 

Lon  °W 

Lat  °N 

Lon  °W 

Length,  km 

AUH 

CMT 

1 

12.5 

-104.4 

10.433 

-103.562 

284.4 

8 

0 

2 

10.04 

-104.33 

8.45 

-104.208 

177 

49 

0 

3 

8.3 

-102.88 

2.667 

-102.15 

630.4 

175 

0 

4 

1.167 

-102.15 

-3.7 

-102.6 

542.5 

142 

2 

5 

-4.088 

-104.522 

-4.481 

-104.644 

45.6 

18 

0 

6 

-4.565 

-106.323 

-6.12 

-106.85 

173 

151 

0 

7 

-6.28 

-107.49 

-8.6 

-108.12 

266.7 

59 

0 

8 

-8.95 

-109.775 

-13.4 

-111.35 

522.9 

21 

0 

Table  2.  Transform  fault  characteristics 


Transform 

Length 

L 

km 

Full  slip- 
rate 

V* 

mm/yr 

Brittle 

area 

A” 

km2 

No. 

AUH 

events 

No.  AUH 
events/ 

A 

km'2 

No.  AUH 
events/ 

L 

km'1 

No.  Of 
Segments 

Segment 
brittle  area 

km2 

Clipperton 

90 

105.3 

326 

129 

0.40 

1.43 

1 

326 

Siqueiros 

140 

111.9 

709.5 

1356 

1.91 

9.69 

5 

277.2 

Quebrada 

120 

137.4 

450.75 

1024 

2.27 

8.53 

4 

163.3 

Discovery 

70 

137.9 

196.7 

586 

2.98 

8.37 

2 

121.4 

Gofar 

190 

138.8 

905.2 

971 

1.07 

5.11 

3 

460.9 

Yaquina 

60 

141.4 

152 

268 

1.76 

4.47 

2 

51.3 

Wilkes 

180 

145 

555 

1505 

2.71 

8.62 

3 

364.9 

Garrett 

130 

149.8 

488 

171 

0.35 

1.32 

4 

179.6 

'DeMets  [1990] 

The  total  brittle  area  of  the  transform  fault  is  the  area  above  the  600°C  isotherm  determined  by  a  half¬ 
space  cooling  model  for  the  given  total  transform  fault  length  (L)  and  spreading  rate  (V). 

Segmented  brittle  area  is  the  sum  of  the  areas  above  the  600°C  isotherm  determined  by  a  half-space 
cooling  model  for  each  segment  of  a  transform  fault  (As  given  in  Table  2). 
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Table  3.  Transform  fault  segmentation  characteristics 


Segment 

name 

Length 

Ls 

km 

Brittle 

area 

As 

km2 

No.  No.  AUH 

AUH  events/ 

events  As 

km'2 

ITSC 

Litsc 

km 

Siqueiros 

Transform 

Si 

33 

64.5 

216 

3.35 

SA 

18 

S2 

42 

98.6 

418 

4.24 

SB 

8 

S3 

24 

38 

201 

5.29 

SC 

5 

S4 

23 

36.5 

237 

6.49 

SD 

5 

S5 

25 

39.6 

284 

7.17 

Quebrada 

Transform 

Q1 

40 

79 

152 

1.92 

QA 

7 

Q2 

24 

34.4 

231 

6.72 

QB 

9 

Q3 

22 

25.9 

175 

6.76 

QC 

10 

Q4 

20 

24 

466 

19.42 

Discovery 

Transform 

D1 

30 

50 

340 

6.80 

DA 

6 

D2 

37 

71.4 

243 

3.40 

Gofar  Transform 

G1 

49 

112 

358 

3.20 

GA 

11 

G2 

35 

63.4 

226 

3.56 

GB 

20 

G3 

89 

285.5 

380 

1.33 

Yaquina  Transform 

Y1 

OSC 

40 

Y2 

31 

51.3 

228 

4.44 

Wilkes  Transform 

WN 

Nannoplate 

1435 

W1 

66 

34.8 

395 

11.35 

WA 

10 

W2 

70 

236.3 

1110 

4.70 

Garrett  Transform 

G1 

18 

18.3 

61 

3.33 

GA 

17 

G2 

34 

59.6 

144 

2.42 

GB 

15 

G3 

22 

25.9 

157 

6.06 

GC 

15 

G4 

40 

75.8 

106 

1.40 
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Table  4.  OSC  characteristics 


OSC 

Offset 

direction 

Offset 

distance  (km) 

Overlap 
distance  (km) 

Overlap 
area  (km2) 

AUH 

# 

CMT 

# 

AUH# 

/km2 

9.0°N 

right 

9 

30 

270 

0 

0 

0 

7.2°N 

right 

2 

- 

2 

20 

0 

10.00 

5.5°N 

right 

8 

13 

104 

9 

0 

0.09 

4.9°N 

right 

5* 

20 

100 

1 

0 

0.01 

3.4°N 

right 

13* 

9 

117 

69 

0 

0.59 

1.0°N 

left 

5* 

- 

5 

- 

0 

0 

1.2°S 

left 

4* 

6 

24 

8 

0 

0.33 

2.8°S 

left 

27* 

21 

567 

65 

2 

0.11 

5.5°S 

right 

15* 

11 

165 

129 

0 

0.78 

7.2°S 

right 

6* 

- 

6 

2 

0 

0.33 

9.6°S 

right 

5* 

10 

50 

2 

0 

0.04 

10.9°S 

right 

6* 

13 

78 

6 

0 

0.08 

*  Data  from  Lonsdale  [1989] 

Table  5.  Teleseismic  and  hydroacoustic  events 


CMT  Event _  _ Hydrophone  Event 


Region 

Mw 

S  L  (dB) 

Year 

Jday 

Time 

Lat 

Lon 

Time 

Lat 

Lon 

Phones 

Siqueiros 

58 

237  22 

1998 

130 

Northern  EPR 

6  35  51  8  93  -103  56 

6  35  50 

8  45 

-102  98 

12456 

Siqueiros 

58 

238  75 

1998 

131 

1  34  40 

8  82 

-103  30 

1  34  42 

8  38 

-103  18 

12456 

Siqueiros 

55 

234  80 

2000 

225 

14  13  15 

8  87 

-103  76 

14  13  11 

8  30 

-104  07 

23456 

Siqueiros 

5  7 

244  59 

2001 

116 

18  59  56 

8  59 

-103  33 

19  00  00 

8  33 

-103  46 

1234567 

Siqueiros 

5  3 

240  90 

2001 

116 

20  39  02 

8  61 

-103  34 

20  39  01 

8  39 

-103  36 

1234567 

Galapagos 

5  3 

203  52 

2001 

141 

16  43  15 

2  67 

-95  13 

16  43  09 

263 

-95  36 

1234567 

2  5SOSC 

54 

231  84 
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Figure  1.  Location  Map  of  the  equatorial  EPR  showing  the  location  of  the  NOAA/PMEL 
Equatorial  Pacific  hydrophone  array  and  major  tectonic  features.  The  locations  of  the 
seven  hydrophones  are  marked  by  white  stars.  White  circles  indicate  the  location  of 
seismic  events  recorded  by  the  AUH  catalog  from  May  1996  to  November  2001 
Overlapping  spreading  center  (OSC)  locations  are  marked  by  white  arrows  (seismically 
active)  and  black  arrows  (no  seismicity).  Inset:  Location  of  the  study  area.  Dashed  black 
lines  and  numbers  illustrate  the  first-order  segmentation  of  the  EPR  within  the  study  area 
(Table  1). 
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Figure  2.  The  Clipperton  and  Siqueiros  transform  systems.  Plate  boundary  is  delineated 
by  the  white  line.  Insert:  The  segmentation  of  the  Siqueiros  transform  is  shown  as  five 
fault  segments,  SI,  S2,  S3,  S4,  and  S5,  and  four  ITSCs,  SA,  SB,  SC,  SD,  (data  from 
Fornari  el  al.  [1989]).  White  circles  indicate  location  of  AUH  recorded  seismicity. 
Error  bars  are  given  by  the  NOAA/PMEL  source  location  error.  CMT  focal  mechanisms 
are  shown  in  red.  Relocation  to  AUH-recorded  event  is  shown  by  red  stars. 
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Figure  3.  The  Quebrada,  Discovery,  and  Gofar  transform  fault  systems.  Plate  boundary 
is  delineated  by  the  white  line,  and  the  geometry  of  the  segmented  systems  (shown  in  the 
inset)  is  labeled  from  north  to  south  from  the  plate  boundary  data  given  by  Searle  [1983], 
White  circles  indicate  location  of  AUH-recorded  seismicity.  Error  bars  are  given  by  the 
NOAA/PMEL  source  location  error.  CMT  focal  mechanisms  are  shown  in  red. 
Relocation  to  AUH-recorded  event  is  shown  by  red  stars. 
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Figure  4.  The  Yaquina  transform  fault  system.  Plate  boundary  is  delineated  by  the 
whilte  line.  The  hybrid  transform  fault  geometry  (shown  at  the  right  panel)  is  from  the 
plate  boundary  given  by  Lonsdale  [1989],  White  circles  indicate  location  of  AUH- 
recorded  seismicity.  Error  bars  are  given  by  the  NOAA/PMEL  source  location  error. 
CMT  focal  mechanisms  are  shown  in  red.  Relocation  to  AUH-recorded  event  is  shown 
by  red  stars. 
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Figure  5.  The  Wilkes  transform  fault  system.  Plate  boundary  is  delineated  by  the  white 
line,  and  the  geometry  of  the  segmented  systems  (shown  in  the  inset)  is  from  the  plate 
boundary  data  given  by  Goff  et  al.  [1993],  White  circles  indicate  location  of  AUH- 
recorded  seismicity.  Error  bars  are  given  by  the  NOAA/PMEL  source  location  error. 
CMT  focal  mechanisms  are  shown  in  red.  Relocation  to  AUH-recorded  event  is  shown 
by  red  stars.  White  boxes  indicate  the  regions  shown  in  Figures  15  and  16. 
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Figure  6.  The  Garrett  transform  fault  system.  Plate  boundary  is  delineated  by  the  solid 
white  line,  and  the  geometry  of  the  segmented  systems  (shown  in  the  inset)  is  labeled 
from  north  to  south  [Hekenian  et  al.,  1992],  White  circles  indicate  location  of  AUH- 
recorded  seismicity.  Error  bars  are  given  by  the  NOAA/PMEL  source  location  error. 
Events  recorded  in  the  CMT  catalog  with  CMT  source  locations  are  given  by  the  focal 
mechanisms.  Relocation  to  AUH  recorded  event  is  shown  by  red  stars. 
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Figure  7.  Northern  OSCs.  Plate  boundary  is  delineated  by  the  white  solid  line.  In  panel 
(b)  dashed  white  line  indicates  the  estimated  plate  boundary  of  Lonsdale  [1985]  based  on 
off-axis  segmentation.  White  circles  indicate  location  of  AUH-recorded  seismicity,  and 
blue  circles  indicate  source  levels  large  enough  to  be  recorded  in  the  REB  catalog.  Error 
bars  are  given  by  the  NOAA/PMEL  source  location  error.  Black  arrows  indicate  the 
direction  of  ridge  propagation. 
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Figure  8.  Southern  OSCs.  Plate  boundary  is  delineated  by  the  white  solid  line.  White 
circles  indicate  location  of  AUH-recorded  seismicity,  and  blue  circles  indicate  source 
levels  large  enough  to  be  recorded  in  the  REB  catalog.  Error  bars  are  given  by  the 
NOAA/PMEL  source  location  error.  Black  arrows  indicate  the  direction  of  ridge 
propagation.  CMT  focal  mechanisms  are  shown  in  red.  Relocation  to  AUH-recorded 
event  is  shown  by  red  stars. 
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Figure  9.  Temporal  and  spatial  statistics  of  the  T-phase  catalog  of  AUH  events,  (a) 
Inter-event  times  of  AUH  events  within  our  study  region:  Small  inter-event  times  indicate 
earthquakes  occurring  closely  spaced  in  time.  Large  inter-event  times  indicate  a  slower 
seismicity  rate.  The  coefficient  of  variation  (Cv)  describes  the  tendency  for  earthquakes 
in  our  study  area  to  cluster  temporally.  Cv  is  defined  as  the  standard  deviation  of  the 
inter-event  times  divided  by  the  mean  inter-event  time.  Values  for  Cv  are  0  for  a  periodic 
data  set,  1  for  Poissonian  (random),  and  >  1  for  clustered  processes  [Kegan  ami  Jackson, 
1991],  (b)  Inter-event  distances  of  AUH  events  within  our  study  region.  Inter-event 
distances  are  grouped  in  bins  of  100  km.  Approximately  55%  of  recorded  earthquakes 
occur  within  0-100  km  of  each  other,  (inset)  Inter-event  distances  grouped  in  bins  of  10 
km.  This  graph  is  zoomed  in  to  the  events  occurring  within  distances  of  250  km  of  each 
other  to  illustrate  the  large  portion  of  events  (nearly  1/3)  that  occur  within  10  km.  Cv  is 
calculated  to  show  the  tendency  of  earthquakes  to  cluster  spatially  within  our  study  area. 
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Figure  10.  The  relationship  between  the  brittle  area  of  the  transform  fault,  given  by  the 
segmented  brittle  area  in  Table  2,  with  (a)  total  seismic  moment  released,  M0,  as  recorded 
in  the  CTM  catalog  and  (b)  total  number  of  recorded  AUH  events.  Linear  fit  is  shown  by 
solid  lines,  with  the  corresponding  R:  fit  values  given. 
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Figure  11.  Number  of  recorded  AUH  events  vs.  transform  fault  segment  length  (Table 

3)- 
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Figure  12.  Earthquakes  recorded  in  the  CMT  for  each  transform  fault  from  1976-2006. 
Event  size  is  scaled  based  on  CMT  recorded  magnitude  (Mw). 
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Figure  13.  Recorded  AUH  events  for  each  transform  fault  over  the  study  period  (1996- 
2002)  plotted  as  their  projected  along-transform  location.  Earthquakes  are  projected  back 
to  the  transform  fault  using  Eulor  pole  location  and  great  circle  rotation  for  each 
transform  fault.  ITSC  location  for  each  transform  fault  is  indicated  by  the  gray  dashed 
lines.  CMT  recorded  events  are  given  by  red  circles,  scaled  to  their  recorded  magnitude. 
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Figure  14.  Time  distribution  of  AUH  events  on  each  of  the  transform  faults  shown  by  1 
month  bins.  (Still  nee  to  add  CMT  events  and  increase  y-axis) 
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Figure  15.  OSC  latitude  vs.  the  timing  of  all  recorded  AUH  events  at  each  OSC  to  show 
rate  and  migration  of  seismicity.  Black  arrows  indicate  the  propagation  direction  for  each 
OSC.  Insets  in  (c)  and  (f)  are  blowups  of  particular  event  sequences  x  and  y, 
respectively,  illustrating  migration  of  earthquakes. 
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Figure  16.  Two  earthquake  sequences  recorded  on  the  Wilkes  nannoplate  boundary:  (b) 
and  (c),  December  6-10,  1998;  and  (d)  and  (e).  May  17-26,  1999.  Plate  boundary  is 
indicated  by  the  solid  white  line  (a,  b,  and  d)  and  solid  black  line  (c  and  e).  Dashed  white 
line  (a,  b,  and  d)  and  dashed  black  line  (c  and  e)  marks  the  location  of  the  Wilkes 
nannoplate  as  discussed  in  Goff  et  ah  [1993],  Earthquake  locations  are  given  by  dots  in 
panels  (b)  and  (d)  with  location  errors  shown.  Panels  (c)  and  (e)  illustrate  the  time 
progression  of  each  earthquake  sequence. 
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Figure  17.  The  June  1-5,  20001  seismicity  sequence  along  the  W2  segment  of  the 
Wilkes  transform  fault  system.  Plate  boundary  is  indicated  by  the  solid  white  line  (a  and 
b)  and  solid  black  line  (c).  (b)  Earthquake  locations,  dots,  with  location  errors  shown, 
(c)  Time  progression  of  the  earthquake  sequence. 
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